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ABSTRACT 

Infrared spectroscopy of the Ha emission lines of a sub-sample of 19 high-redshift 
(0.8 < z < 2.3) Molonglo quasars, selected at 408 MHz, is presented. These emission 
lines are fitted with composite models of broad and narrow emission, which include 
combinations of classical broad-line regions of fast-moving gas clouds lying outside the 
r ) ' quasar nucleus, and/or a theoretical model of emission from an optically-thick, flat- 

tened, rotating accretion disk, with velocity shifts allowed between the components. 

■ All bar one of the nineteen sources are found to have emission consistent with the 

presence of an optically-emitting accretion disk, with the exception appearing to dis- 

q , play complex emission including at least three broad components. Ten of the quasars 

!h ' have strong Bayesian evidence for broad-line emission arising from an accretion disk 

r/y \ together with a standard broad-line region, selected in preference to a model with 

& . two simple broad lines. Thus the best explanation for the complexity required to fit 

the broad H a lines in this sample is optical emission from an accretion disk in ad- 
dition a region of fast-moving clouds. We derive estimates of the angle between the 
rotation axis of the accretion disk and the line of sight. Deprojecting radio sources 
■^j- ' on the assumption of jets emerging perpendicular to the accretion disk gives rough 

f^) , agreement with expectations of radio source models. The distribution in disk angles 

■ is broadly consistent with models in which a Doppler boosted core contributes to the 

chances of observing a source at low inclination to the line of sight, and in which the 
radio jets expand at constant speed up to a size of ~ 1 Mpc. A weak correlation is 
found between the accretion disk angle and the logarithm of the low-frequency radio 

^\ \ luminosity. T his is direct, alb eit tenuous, evidence for the receding torus model first 

C^) . suggested by iLawrencd ([T991) in which the opening angle of the torus widens with 

increasing radio luminosity. The highest accretion disk angle measured is 48°, consis- 
tent with the opening angle predicted for radio-luminous sources. Velocity shifts of 
the broad Ha components are analysed and the results found to be consistent with 
a two-component model comprising one single-peaked broad line emitted at the same 
redshift as the narrow lines, and emission from an accretion disk which appears to be 
preferentially redshifted with respect to the narrow lines for high-redshift sources and 
blueshifted relative to the narrow lines for low-redshift sources. An additional analy- 
sis is performed in which the disk emission is fixed at the redshift of the narrow-line 
region; although only two quasars show a robust change in fitted angle, the radio lu- 
minosity - disk angle correlation falls sharply in probability, and so is strongly model 
dependent in this sample. 

Key words: galaxies: active - galaxies: jets - quasars: emission lines - quasars: 
general. 
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1 INTRODUCTION 
1.1 Orientation effects 

Optical emission spectra of Active Galactic Nuclei (AGN) 
matched in radio and optical luminosity are now firmly be- 
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lieved to be strongly influenced by orientation effects, with 
only small underlying differences in the sources themselves. 
There are two separate orientation-dependent effects which 
alter the optical spectra of AGN. 

The Type 1/Type 2 classification of AGN is made ac- 
cording to the presence of broad emission lines. Type 2 
AG N possess only narrow emission lines, < 2000 km s _1 
(e.g. |Petersonlll997l ) and weak non-stellar continuum emis- 
sion. Type 1 AGN have broad emission lines of ~ 2000 
- 20000 km s , and strong non-stellar continuum emis- 
sion, in addition to narrow lines similar to the Type 2 
sources. An explanation for this disparity grew from the 
discoveries of broad emission li nes seen in polarised ligh t 
from Type 2 Seyfert sources fe.g. lAntonucci fc Millerll 19851 ). 
which suggested orientation-dependent obscuration caused 
by an inter vening screen of matter, s uch as a dusty molec- 
ular torus l|Krolik fc Begelmanl Il986l ). It has recently be- 
come clear that obscuration by dust in starbursting galaxies 
can also be responsible for c oncealing Type 2 AGN (e.g. 
Mar tmez-Sansigre et al.ll2005l ) . 

The second orientation-dependent effect arises from 
the relativistic moti o n of the plasma in the radio jets. 
IScheuer fc Readheadl (|l979l ) first suggested that viewing a 
radio source with opposing relativistic jets would cause a 
large contrast in the luminosities of the approaching and re- 
ceding jets, and that objects with jet axes close to the line 
of sight would b e seen more of t en du e to Doppler boost- 
ing of the core. lOrr fc Browne! l|l982l ) made the connec- 
tion between Doppler boosting of the core, and a measure 
of quasar orientation from the core-to-lobe radio flux den- 
sity ratio; this allowed them to unify the "core-dominated" 
quasars with flat optical spectra, viewed at angles close to 
the line of sight, with the "steep-spectru m" or radio-lobe- 
domin ated quasars viewed at larger angles. Wills fc Brownd 
|l986t ) linked the radio properties to the optical properties 
by discovering an anticorrelation between the core-to-lobe 
radio flux ratio and the width of the H/3 line, interpreting 
this connection as the result of beaming of radio emission 
from a jet emerging along the rotation axis of an accretion 
disk; the broad H /3 lines arise from the accretion disk with 
a width correlated with the angle between the line of sight 
and the disk axis. 

The two optical schemes were reconciled by iBarthell 

(1989), who gave a consistent picture in which FRII nar- 
row line radio galaxies, steep-spectrum quasars and flat- 
spectrum quasars are all drawn from the same parent pop- 
ulation, but viewed at decreasing angles to the line of sight. 
A revie w of these so-called unifie d s chemes for AGN c an be 
found in lUrrv fc Padovanil l|l995l ) or I Antonucc] j 19931 ). 

A solid understanding of how quasar emission lines 
arise, and how they are affected by the AGN environment 
along different sight lines, is not only an interesting study in 
terms of quasar structure, but is also vital in order to dis- 
entangle orientation effects from large-scale AGN surveys 
which enable the study of cosmic evolution. 

1.2 Accretion disks 

There is a growing body of evidence that AGN are pow- 
ered by accretion of gas and dust onto supermassive black 
holes, and this is now the accepted paradigm. As the 
black hole feeds on the surrounding material, it is expected 



that this will form an ac cretion disk of infalling matter 
l|Shakura fc SunvaevllT973T ). The current theory is that ac- 
cretion disks have two parts: a puffcd-up, X-ray-emitting 
inner disk, and a flattened, outer disk which gives rise to 
bro ad optical emission lines. 

ICollin-Souffrin et al. I l|l980l ) first suggested that a thick- 
ened inner accretion disk, shielding and reprocessing the 
hard X-ray emission from the black hole, gives rise to 
the optical low-ionisation Fell emission seen in Type 1 
Seyferts, from the atm osphere above the geometrically-thin 
outer part of the disk. iRees et"all l|l982l ) postulated that a 
geometrically- and optically-thick ion-supported torus sur- 
rounds the supermassive black hole at the ce ntre of a radio 
galaxy , collimating the emerging radio jets. iTanaka et alj 
( 1995) observed a broad, asymmetric iron Ka line consis- 
tent with emission from a disk of this description, situated 
between ~ 3 - 10i?G (where _Rg = GM/c 2 is the gravita- 
tion al radius) from th e nucleus of the AGN. 

iFilippenkol (| 198ST ) reviewed the different lines of evi- 
dence from optical and UV data that flattened, extended 
accretion disks fuel the central black holes of galaxies. For 
example, iBaldwinl (|l977l ) recorded the anticorrelation of the 
UV continuum luminosity with the equivalent widths of 
broad Civ emission (the "Baldwin Effect"). Both this ob- 
servation and the "bi g blue bump" of excess UV contin- 
uum emission found bv lMalkan fc Sargentl (|l982l ) may be ex- 
plained by emission from an opt ically-thick, geometrically- 
thin accretion disk (jNetzerlll985l ). The strongest direct ev- 
idence for this scenario is the presence of double-peaked, 
low-ionisation optical emission lines seen in some AGN, 
which arise from the Doppler effect acting on the emis- 
sion from rotating mat er ial in the ou t er acc retion disk (e.g. 
IChen fc Halpernl ll 19891 ). iPerez et all (|l988l )). 

The thin, optically-emitting disk must be illuminated 
by some mechanism. The photoionising flux may originate 
either from the central no n-thermal source, or from t he in- 
ner, X-ray-emitting disk (|Collin-Souffrin et al.| [l980): and 
the radiation may illuminate the outer disk directly, or be 
scattered f rom a highly-ioni sed diffuse medium above the 
outer disk (|Chen etal.ll 19891 ). 

Optical double-peaked lines have to date only been 
found in a r elatively low percent a ge of ra dio-loud AGN 
(~ 10 %, see lEracleous fc Hafpernl l)l994 )). IStrateva et~aD 
(|2003h discovered that radio-quiet quasars also emit double- 
peaked lines, although these appear to be rarer still: double- 
peaked emission was seen in ~ 3% of the SDSS AGN, includ- 
ing both radio-loud and radio-quiet s ources. Double-peaked 
profiles are not unique to Balmer lines: IStrateva et al.l (|2003l ) 
found double-peaked Mgn emission lines in some SDSS 
AGN. 

It is not clear why the double-peaked line profiles should 
be rare, although there are several possibilities: the outer ac- 
cretion disk may simply be obscured by broad-line-emitting 
clouds surrounding it; or if the outer accretion disk causes 
a wind, then the broad lines w hich arise from this wind are 
predicted to be single-peaked ()Murrav fc Ch iang 1993). In 
any case, it should not be expected that the low-ionisation 
broad lines seen in an AGN originate solely from the rotat- 
ing disk; single-peaked emission may be seen in addition to 
double-peaked profiles. 

The accretion dis k mod el used in this paper is taken 
from IChen fc Halpernl §989), and consists of a thick, hot 



© 2008 RAS, MNRAS QOO.[Tti36l 



Modelling the orientation of accretion disks in quasars using Ha emission 3 



torus, whose outer edge may reach up to 100 Rg from the 
black hole. Inverse Compton-scattered X-rays from this in- 
ner disk illuminate an optically-thick, flattened outer disk 
l|Halpern fc Chenlll989l ). The thin, circular disk, which pro- 
duces the double-peaked emission lines, may extend up to 
~ 10 s Rg from the central engine. The distinctive line pro- 
files are caused by the rotation of the disk, splitting the 
emission into redshifted receding material and blueshifted 
approaching material; the blueshifted peak is of higher inten- 
sity than the redshifted peak, as a result of Doppler boost- 
ing. 

The chosen model was necessarily simple, to limit 
the number of free parameters. More complex disks, such 
as an elliptical disk, which may ar ise when a single star 
is di srupted near the black hole (jGurzadvan fc Ozernovl 
Il979l ). or a wa rped disk, tho ught to occur around rotat- 
ing black holes (|Bachevlll999h . give rise to a wider range of 
line profiles, including double-peaked profiles with a red- 
ward peak of highe r intensity than the blueward peak. 
IStrateva et al.l (2003) found that assuming all low-ionisation 
broad-line emission comes from a disk, non-axisymmetric 
disks would be required in ~ 60 % of their SDSS sample, 
while lEracleous fc Halpernl (120031 ) found that in their sam- 
ple of 106 radio-loud AGN, 20% have double-peaked Ha 
profiles visible to the eye, of which ~ 40% require a model 
more complex than the circular Keplerian disk. 

In this paper, a small, but close to complete, sample of 
radio-loud quasars are analysed to determine if their spec- 
tra include emission from circular, planar accretion disks, 
either as the sole component of broad optical emission, or in 
combination with a single-peaked broad line. 



1.3 Velocity shifts 

It has been generally accepted that the narrow-line region 
(NLR) of an AGN falls near the systemic redshift. The NLR 
is extended, and appears to be reasonably independent of 
viewing angle effects, and so the narrowness of the lines con- 
strain s the velocity of this gas to be small. Irleckman et al.l 
(jl98ll ) showed, for a handful of low-z Seyferts and radio 
galaxies, that the narrow lines had small blueshifts of be- 
tween ~ - 300 km s _1 relat ive to the neutral hydro - 
gen emission of the host galaxies. IVanden Berk et al.l ()200ll ) 
discovered that, for composite spectra created with several 
thousand Sloan Digital Sky Survey (SDSS) quasar spectra 
covering a wide redshift range (0.04 < z < 4.8), the nar- 
row lines are blueshifted by small velocities (< 100 km s _1 ) 
which correlate with ionisation potential. 

The broad-line region (BLR) has a complex struc- 
ture, with emission line redshifts which vary according to 
species, implying separate regions of gas (see Figure 6 of 
ICollin-Souffrin et all (ll980l )). It is thought that the high- 
ionisation lines (HILs) arise from a compact, spherical region 
close to the central black hole, while the low-ionisation lines 
(LILs) are formed in a flattened struct ure further from th e 
AGN centre, possibly an ac cretion disk (jKrolik et al.lll99lT ). 
ICollin-Souffrin et all (|l988r i suggested that the HILs might 
be produced by sho cks in an outflowing wind. 

iGaskelll (|l982l ) demonstrated for a sample of flat- 
spectrum quasars with z ~ 0.2 - 2.3 that the HILs, such as 
Cm], Civ and Nv, are blueshifted by ~ 600 km s -1 with 
respect to the LILs, which include Mg 11, O I and the Balmer 



lines. IWilkesI (|1984| ) confirmed this shift for high- z quasars 
(2 ^ z < 3), finding a slightly higher range of shifts, up to ~ 
1400 km s _1 . Blueshifted HIL zones have also been observed 
bv lEspev et~aH (Il989l ). who found shifts of ~ 1000 km s" 1 
in a s mall sample of 1.3 < z < 2.4 sources, and ICorbinl 
(1990), who found shifts exceeding 4000 k m s -1 for lumi- 
nous, optically-selected sources with z > 1. 1 Richards et all 
(2002) also confirm this trend from studies of ~ 800 quasars 
with 1.5 < z < 2.2 from the SDSS, though they find a wide 
distribution of velocity shifts of the high-ionisation C IV with 
respect to the low-ionisation Mgn, ranging from redshifts of 
~ 500 km s _1 to blueshifts of over 2000 km s~\ and they 
take pains to point out that they do not believe it is a line 
shi ft so much a s a lack of flux in the red wing of the line. 

iMacIntoshl l| 19991 ) found that the HIL zone is at the 
same redshift as the narrow lines in low-z sources, but 
for a sample of quasars with 2 < z < 2.5, broad H/3 
had a redshift of ~ 500 km s" 1 relative to [Oiii]. Al- 
though the LILs are typically considered to have low ve- 
locity shifts, there have been recorded instances of Balmer 
lines with very high redshifts, e.g. ~ 2100 km s" 1 for 3C277 
llOsterbrock. Koski fc Phillipsjll976h and ~ 2600 km s" 1 for 
OQ208 |Osterbrock fc Cohenll 19791) . 



2 SAMPLE SELECTION 

A sub-sample of 1 9 quasars was defin ed from the Molon- 
glo Quasar Sample |Kapahi et al.l [l998) . The Molonglo sam- 
ple of radio sources was selected from the Molonglo Refer- 
ence Catalog (MRC) l|Large et al.lll98ll ). a 408 MHz survey 
conducted with the Molonglo Synthesis Telescope which is 
99.9% complete at 1 Jy. The radio sources were then identi- 
fied from VLA 1-arcsecond resolution r adio images, optica l 
imaging and spectroscopy, the quasars (Kapahi ct al. 1998) 
being distinguished from the radio galaxies ()McCarthv et all 
1 19961 ) by the presence of broad optical emission lines. 

The Molonglo quasar sample includes all quasars with 
flux densities 5408MHz > 0.95 Jy in a 10°-wide strip in the 
Southern sky, —20° > S > —30°, excluding sources with low 
Galactic latitu de (| b |> 20°) and a lso a strip in Galactic 
R.A. (details in lKapahi et alJ (|l998l )) to define a sample of 
manageable size. It should be noted that as the Molonglo 
sample was selected at the mid-range frequency of 408 MHz, 
there are likely to be some sources included in the sample by 
virtue of their strong radio core emission, and are therefore 
inclined at small angles to the line of sight; this builds an 
orientation bias into this survey. 

The quasar sub-sample was selected on the basis of four 
criteria: observability during the relevant time period; red- 
shift such that H a and H /3 emission falls in wavelength 
windows corresponding to regions of high atmospheric trans- 
parency (0.8 < z < 1.0, 1.5 < z < 1.65, 2.2 < 2 < 2.3); 
sufficient J-band brightness to be observable in a reasonable 
integration time; and exclusion of the RA range 14h - 03h 
in accordance with scheduling constraints. The J-band mag- 
nitude limit is the only factor which adds a significant bias 
in the sub-sample selection. The objects were selected to be 
brighter than J ~ 18.5, and this excluded one source from 
the sample, MRC0418-288. This source is likely to be red- 
dened or dusty, which means that it has a higher chance of 
being inclined at a large angle to the line of sight. MRC1256- 
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243 is an extra core-dominated source which was not ob- 
served due to scheduling limitations; this source is likely to 
have been boosted into the sample by virtue of its strong 
core emission in any case. 

Throughout the paper, a cosmology of Hq = 70 km s _1 
Mpc -1 , Qm = 0.3 and Qa = 0.7 is assumed. 



3 INFRARED SPECTRA 
3.1 Data acquisition 

Near infrared spectra of H a were obtained for the Molonglo 
quasar sub-sample with the Infrared Spectrometer And Ar - 
ray Camera (ISAAC) spectrograph (|Moorwood et allll998l ) 
at ESO's VLT UT1, in service mode, between October 
2001 and February 2002. The spectra were taken in short- 
wavelength, low-resolution mode, using a 1-arcsec-wide slit. 
The observational details, including exposure times, seeing 
and airmass, are given in Table [T] 



3.2 Data reduction 

The raw images were cleaned of cosmic rays in two stages: 
crmedian in IRA$\ was used to automatically remove cos- 
mic rays which fell more than 10cr below or 3<r above the 
median value, replacing these pixels with the median value. 
The regions of the image including the spectrum and the sky 
subtraction zone were then cleaned by hand using credit in 
IRAF, replaci ng bad pixels with local sky values. Known de- 
tector effects (|Amico et alj r20021 were corrected as follows. 
"Electrical ghosts" are additional signals i n the image, and 
were removed using the dedicated eclipse l|Devillardlll997n 
recipe ghost. The "odd-even column effect" causes an inten- 
sity difference in the alternate rows of the image; this was 
cleaned from images in which it was apparent using a script 
which masked the two pixels corresponding to the highest 
frequency variations in Fourier sp ace for each quadra nt of 
the image, following the method of I Amico et alJ (|2002l ). 

The images were flatfielded, corrected for distortion ef- 
fects and wavelength calibrated using arc lamp frames for 
the relative calibration and OH sky emission lines to zero- 
point correct the calibration; these procedures were carried 
out using standard techniques in IRAF. The sky background 
is strong and time-variable in the infrared, and so ISAAC 
spectra are observed in nod-and-jitter mode, producing pairs 
of images with the spectra located on different regions of the 
CCD. For each image, the sky background was subtracted 
using the neighbouring image, before all the spectra were 
stacked using a script in IDL. 

The spectra were extracted with the IRAF routine 
apall, and telluric features were removed and a flux calibra- 
tion applied simultaneously using one reference star for each 
quasar, with IDL procedures. Finally, the spectra were cor- 
rected for the dust reddening of the Milky Way, using the 



1 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



ISchlegel et all (|l998T l maps of Galactic dust emission, and 
the spectra were corrected to the heliocentric rest-frame 0- 

3.3 Molonglo infrared spectra 

Spectra of the quasars are presented in Figure [1] For the 
purposes of this figure, these have been smoothed in D/PSClfl 
with a Gaussian filter of FWHM 4.7 pixels (pixel scales vary 
slightly between spectra, but are typically 2 A pixel -1 ) to 
reduce the noise. 

3.4 Notes on infrared spectra 

MRC0222-224: This spectrum contains noise spikes from 
poor sky subtraction which affects most of the H-band spec- 
tra. The apparent absorption line in the broad Ha line is 
an artifact, and is excluded from the emission line fitting. 
MRC0327-241: The continuum slopes with an index of 
a ~ 1.7, where F v oc v~ a . This is a BL-Lac typ e continuum, 
as the se have spectral indices of a > 1 (e.g. iBrown et al.l 

dUi)). 

MRC0346-279: The continuum is strongly sloped, with 
a ~ 2.7; this is a BL-Lac type continuum. 
MRC0413-210: The narrow line at a wavelength of 7322 
A (observed 13230 A) is O II. This line was not fitted in the 
analysis, as it is well-separated from the Ha region. 
MRC0413-296: The apparent narrow features to the right 
of the broad H a profile are artifacts of the sky subtraction. 
MRC0430-278: This spectrum has a large number of noise 
spikes remaining from the sky subtraction process. 
MRC0450-221: The drop in flux longwards of 1350 A is 
not a real feature, but an artifact introduced during the flux 
calibration. This part of the spectrum is excluded from the 
emission line fitting. 

MRC1019-227: The sky subtraction is poor, leading to 
artifacts which resemble narrow lines. The apparent narrow 
line at the wavelength of Ha is a sky line defect, and is 
excluded from the emission line fitting. 

MRC1217-209: This spectrum is very noisy due to poor 
seeing. The structure of Ha is not readily apparent. 
MRC1301-251: Broad emission is weak in this spectrum 
compared to the strong narrow lines. 

MRC1355-215: The apparent narrow absorption features 
are artifacts from imperfect sky subtraction. 

3.5 Redshift measurements 

Improved redshifts for these quasars were measured from 
the ISAAC spectra, or from new optical spectra (Janssens 
et al., in preparation) in cases that strong narrow lines were 
available ([Oiii] for preference, followed by Balmer lines). If 
no strong narrow lines were present in the new spectra, then 
the most recent measurements from the literature were se- 
lected. The only exception to this was MRC1349 - 265, w hose 
redshift was given as z = 0.934 in lBaker et"al] <|l999h . but 
for which z = 0.924 was a solid measurement, even in the 

2 The correction was calculated using ephemerides from the 
Markwardt IDL library. 

3 DIPSO is a Starlink program 
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Figure 1. Observed near-infrared spectra, including the Ha emission lines, for the Molonglo sub-sample quasars, observed 
with ISAAC on ESO's VLT UT1. The observed wavelength is in A, and the flux density scale is in mjy. The spectra have 
been slightly smoothed. 
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3 


MRC0346-279 


03 


18 


38.1 


-27 


49 


14 


0.989 


20.5 


2001-10-12 


J 


180 


x 8 


0.7 


1 


077 


i 


MRC0413-210 


04 


16 


04.3 


-20 


56 


28 


0.807 


18.4 


2001-10-12 


J 


180 


x 8 


0.9 


1 


053 


5 


MRC0413-296 


04 


15 


08.7 


-29 


29 


03 


1.614 


18.6 


2001-10-12 


H 


180 


x 8 


0.7 


1 


004 


6 


MRC0430-278 


04 


32 


17.7 


-27 


46 


24 


1.633 


21.3 


2001-12-23 


H 


180 


x 10 


1.2 


1 


138 


7 


MRC0437-244 


04 


39 


09.2 


-2d 


22 


08 


0.834 


17.5 


2001-10-06 


J 


180 


x 8 


0.6 


1 


012 


8 


MRC0450-221 


04 


52 


44.7 


-22 


01 


19 


0.900 


17.8 


2001-11-20 


J 


180 


x 8 


0.5 


1 


214 


9 


MRC0549-213 


05 


51 


58.3 


-21 


19 


49 


2.245 


19.1 


2001-12-22 


K 


100 


x 28 


0.5 


1 


075 


10 


MRC1019-227 


10 


21 


27.6 


-23 


01 


54 


1.542 


21.1 


2001-12-25 


H 


180 


x 10 


0.7 


1 


144 


11 


MRC1114-220 


11 


16 


54.5 


-22 


16 


53 


2.286 


20.2 


2002-01-01 


K 


100 


x 28 


0.7 


1 


137 


12 


MRC1208-277 


12 


10 


43.6 


-27 


58 


55 


0.828 


18.8 


2002-01-27 


J 


180 


x 8 


0.5 


1 


034 


13 


MRC1217-209 


12 


20 


22.3 


-21 


13 


09 


0.814 


20.2 


2002-02-05 


.1 


180 


x 8 


1.1 


1 


106 


14 


MRC1222-293 


12 


25 


01.2 


-29 


38 


17 


0.816 


18.5 


2002-01-27 


.1 


180 


x 8 


0.6 


1 


006 


15 


MRC1301-251 


13 


04 


14.7 


-25 


21 


37 


0.952 


21.0 


2002-02-12 


J 


180 


x 8 


0.1 


1 


Oil 


16 


MRC1349-265 


13 


52 


10.3 


-26 


49 


28 


0.924 


18.4 


2002-02-11 


J 


180 


x 8 


0.5 


1 


051 


17 


MRC1355-215 


13 


58 


38.2 


-21 


-18 


54 


1.607 


19.9 


2002-02-13 


H 


180 


x 8 


0.5 


1 


002 


IS 


MRC1355-236 


13 


58 


32.7 


-23 


52 


20 


0.832 


17.8 


2002-02-11 


J 


180 


x 8 


0.4 


1 


013 


19 


MRC1359-281 


11 


02 


02.4 


-28 


22 


25 


0.802 


18.7 


2002-02-14 


.1 


180 


x 16 


0.6 


1 


048 



Table 1. Observational details: Columns 1 and 2: Quasar index and Molonglo Reference Catalog (MRC) name; Columns 3 
and 4- J(2000) right ascension and declinat ion; Column 5: Redsh ift; Column 6: Optical bj magnitude from the UK Schmidt 
IlaJ plates, where bj = B - 0.23(B - V) l|Bahcall fc boneir al ll98Gv i: Column 7: Date of observation; Column 8: Waveband of 
observation; Column 9: Exposure time in seconds multiplied by the total number of exposures in the observation; Column 
10: Seeing in arcseconds, measured from the images; Column 11: Airmass. 

at m 



absence of narrow [O ill]; a change of this size is likely to re- 
sult from a typographical error in the original paper. Table 
[2] lists the redshifts and their origins, in addition to mea- 
sured FWHMs of the broad lines, the integrated fluxes of 
the broad and narrow H a lines, and core-to-lobe flux ratios 
at 10 GHz in the rest frame (-Rioghz). 



4 MODELS OF EMISSION 
4.1 Set of models 

The emission in the rest-frame range ~ 6000 - 7000 A was 
modelled as the sum of narrow line emission and broad emis- 
sion. There are four sets of models with different broad emis- 
sion contributions. One set of models has a single Lorentzian 
line representing emission from a classical BLR of fast- 
moving clouds; one set comprises a Lorentzian and a Gaus- 
sian line, to simulate two BLRs at different temperatures, 
or a BLR plus an outflow. There are two sets of models in- 
cluding accretion disks: one set has the accretion disk plus 
a Lorentzian profile to represent a standard BLR, while the 
other set includes broad emission from the accretion disk 
only. There are three models in each set: one with narrow 
H a plus O i, S n and N n lines; one with narrow H a but 
none of the forbidden lines; and one with no narrow line 
emission. There are therefore twelve models constructed in 
this modular way, as shown in Figure [5] The parameters in- 
cluded in each of the models are detailed in Table [3] and the 
prior probability ranges for these parameters are shown in 
Table H 




o ra — 

-z- -z- < 



Lorentzian plus disk emission 


5 


3 


1 


Lorentzian plus Gaussian line 


9 


8 


7 


Disk emission only 


12 


11 


10 


Single Lorentzian 


6 


4 


2 



Figure 2. An illustration of the components included in each of 
the models, 1 to 12. 



4.2 Continuum emission 

Broadband emission from the central engine contributes a 
smooth continuum to the spectrum; this part of the emis- 
sion was not included in the models, since it is dominated by 
the processes occurring immediately around the black hole, 
and not the distribution of gas and dust outside the central 
region. Instead, it was subtracted from the spectrum, using 
either a linear or a quadratic fit made in DIPSO. The ne- 
cessity of the quadratic term was judged by eye. In order to 
allow for a small residual component of continuum emission, 
a constant term was included in each of the models. 
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Quasar 


z 


Source 


FWHM 


^Narrow H a 


^Broad H a 


Measured 


RiOGHz from 


Source of 








of z 


(km s" 1 ) 


(10- 19 W m~ 2 ) 


(10- 19 W m~ 2 ) 


^lOGHz 


literature 


^?10GHz 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


1 


MRC0222-224 


1.603 


2 


2230 


46 


740 


0.0065 


< 0.37 


1 


2 


MRC0327-241 


0.895 


3 


3420 


5 


255 


10.0 


> 1 


1 


3 


MRC0346-279 


0.989 


2.3 


2320 


ND 


144 


10.2 


> 5 


1 


■1 


MRC0413-210 


0.807 


2.3 


2250 


58 


1054 


0.78 


0.71 


1 


5 


MRC0413-296 


1.614 


1 


5530 


71 


2573 


0.020 


0.031 


1 


6 


MRC0430-278 


1.633 


1 


2480 


4 


581 


0.46 


— 




7 


MRC0437-244 


0.834 


2 


4950 


30 


3477 


0.10 


0.098 


1 


8 


MRC0450-221 


0.900 


2 


6620 


57 


3739 


0.060 


0.086 


1 


9 


MRC0549-213 


2.245 


2,3 


2630 


ND 


561 


0.31 


< 0.55 


1 


10 


MRC1019-227 


1.542 


1 


1910 


6 


377 


0.051 


— 




11 


MRC1114-220 


2.286 


1 


2650 


55 


1715 


0.078 


0.08 


2 


12 


MRC1208-277 


0.828 


3 


3490 


144 


1886 


0.049 


0.088 


1 


13 


MRC1217-209 


0.814 


2.3 


5850 


ND 


537 


0.078 


0.057 


1 


14 


MRC1222-293 


0.816 


3 


2840 


77 


935 


0.98 


0.38 


1 


15 


MRC1301-251 


0.952 


3 


1530 


52 


343 


0.014 


0.020 


1 


1(5 


MRC1349-265 


0.924 


1 


3350 


52 


2366 


0.25 






17 


MRC1355-215 


1.607 


1 


2600 


25 


563 


0.36 


0.84 


1 


18 


MRC1355-236 


0.832 


3 


2180 


76 


1007 


0.097 


0.098 


1 


19 


MRC1359-281 


0.802 


2.3 


3820 


25 


575 


0.12 







Table 2. Measured properties of the quasar sub-sample: Columns 1 and 2: Quasar index and MRC name: Colum n 3 : Redshift; Col- 
umn 4 : Origin o f redshift measurement - 1 = New measurement presented in this paper, 2 = lBaker et"ai] fll999h . 3 = lde Silva et al.l 
(in preparation); Column 5: FWHM in km s~ 1 ; Column 6: Narrow Ha line flux in units of 10 — 19 W m — 2 . ND in this column 
indicates that narrow Ha was not detected. Column 7: Broad Ha line flux in units of 10~ 19 W m~ 2 . Column 8: Measured 
core-to-lobe flux ratio at rest-frame 10 GHz (i?ioGHz); Column 9: Core-to-lobe flux rat io at rest-frame 10 GHz ( -RioGHz) from 
the literature; C olumn 10: Origin of the core-to-lobe flux ratio from the literature — 1 = iKapahi et al. (1998), 2 = Ide Silva et alj 
(|in preparation! ). 



4.3 Narrow-line emission 

All narrow lines were modelled with Gaussian distributions. 
The relative intensities of some of the forbidden lines were 
constrained, either because they are fixed by transition prob- 
abilities, or because they depend on temperature and elec- 
tron density, which can be assumed to be approximately 
constant within the NLR. The Nil and Ol emission lines 
are temperature sensitive, and may depend on other factors 
such as reddening; however, since these emission lines make 
only a modest contribution to the spectra, the line ratios of 
O I6300 /0 16364 and N H6583 /N H654 8 were fixed at a value of 
3.0, with reference to iKoskil (| 19781 ). The S Ilene/S Her3i ra- 
tio depends both on the square root of the temperature and 
on the electron density; this ratio was there fore left as a fre e 
parameter, allowed to vary from 0.2 - 3.0 l|Petersonlll997t ). 
Reasonable prior ranges of 0.003 - 10 for the relative inten- 
sities of S1I6732, N1I6550 and 16302 with respec t to na rrow 
H a were estimated from IVeilleux fc Osterbrockl (|l987l ) . 

The widths of the observed narrow lines depend upon 
the effective spectral resolution. The intrinsic line width due 
to Doppler broadening of narrow Ha was left as a free pa- 
rameter in the model, and assumed to be the same for all 
narrow lines, since they are all low-ionisation lines formed 
in approximately the same region. This intrinsic width was 
convolved with the line width due to the spectral resolution, 
which is a wavelength-dependent quantity. The spectral res- 
olution was measured from the night sky OH lines for each 
waveband, and was found to be in the range 400 - 600. In 
cases where the quasar did not fill the spectrograph slit, this 
resolution was scaled down by the ratio of the seeing (esti- 



mated for each spectrum by its spatial extent) to the slit 
width. 



4.4 Broad-line emission 

Emission from the BLRs was modelled as a Lorentzian line, 
which is collisionally broadened with a width proportional 
to P/T, where P is the pressure and T the temperature. 
This line profile was chosen since AGN broad lines have been 
found observationally to have broader wings than Gaussians 
(e.g. lPetersoiil (|l997f )). and Lorentzian profiles fulfil this re- 
quirement, whilst being simple and smooth. 

One set of the models has an additional broad line com- 
ponent, which can represent a range of different physical 
processes, including two separate BLRs with different tem- 
peratures, or outflows from the central region. A red wing 
on the emission line may be caused by an outflow of op- 
tically thick clouds, in which case only those on the far 
side of the quasar f rom the central source w ould be visi- 
ble l|Capriotti et al.1 1| 19791 ). ISmith et all (|l98ll )). This broad 
component is modelled by a Gaussian profile, which arises 
from Doppler or thermal broadening, and so has a width 
proportional to VT. The difference between Lorentzian and 
Gaussian profiles is minimal near the line centroid, though 
the Lorentzian has much broader wings; combining the two 
different lines allowed the maximum degree of flexibility in 
the two-component BLR models. 
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Parameter 

1 


2 


3 


4 5 
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Broad H a Lorentzian width • 
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Broad H a intensity • 
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Narrow H q central wavelength • 
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# 


• 
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7 


Narrow H q Gaussian width • 


• 
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# 


• 




• 


• 


8 


Narrow H a intensity • 




• 
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• 


• 
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N H5550 /Narrow H ct intensity ratio • 














• 




10 


S Il67i8 /S 115732 intensity ratio • 














# 




11 


S H6732 /Narrow H o. intensity ratio • 


















12 


16302 /Narrow H ol intensity ratio • 








# 






• 




13 


Second broad Ha component shift with re- 










• 


• 








spect to narrow Ha 


















14 


Second broad H a component Gaussian width 










• 


• 






15 


Second broad H a component intensity 










• 


• 







16 Disk intensity normalisation • • • • • • 

17 Disk shift with respect to 6564.61 A • • • • • • 

18 Sine of the disk angle • • • • • • 

19 Local velocity dispersion of the disk material • • • • • • 

20 Inner disk radius • • • • • • 

21 Multiplication factor for disk outer radius • • • • • • 

Number of free parameters 17 11 13 7 10 4 14 10 7 14 10 7 



Table 3. Summary of the components included in each model. Columns 1 and 2: Index and description of the free 
parameters; Columns 3 - 12: points indicating that a parameter is included in a particular model. Figure[2]gives a simplified 
summary of this information. 



4.5 Accretion disk emission 

The template used for the ac cretion disk emission was taken 
from lChen fc Halpernl ij 19891 ). and describes emission from a 
geometrically-thin, optically-thick disk, which is illuminated 
by a thick, hot inner disk. The model assumes that the disk 
is circular and Keplerian. In a disk such as this, viewed at a 
non-zero angle to the line of sight, the receding material is 
redshifted, and the approaching material blueshifted; addi- 
tionally, the blue peak has a higher intensity due to Doppler 
boosting. If the accretion disk is viewed face on (at an in- 
clination angle of zero) , there is no velocity difference along 
the line of sight, and the disk emission is single-peaked. As 
the angle increases, the velocity of the disk material along 
the line of sight increases, and the red and blue peaks move 

further apart. 

Following [Ch en fc Halpernl l|l989h . the expression for 
the disk emission per unit frequency interval, F v , is 

F » = K j R J K ex P~ 2D^Av 2 p3fil ~M-P)dfl#(l) 
where 

y G 2 26 M 2 cos6> diak 
and 



ff (D) = l+I/ — j l\ (3) 

fil^cos 2 disk + J R(I?-(l-|)5) 2 I 

and G is the gravitational constant, c is the speed of light; 
M is the black hole mass; eo is the normalisation and q is 
the radial exponent of the disk emissivity, defined by e = 
(eo/47r).R -9 ; d is the luminosity distance of the quasar; #di s k 
is the disk angle, defined as the angle between the accretion 
disk rotation axis and the line of sight; vq is the rest-frame 
frequency of the line emission; Av is the dimensionless local 
velocity dispersion of the disk material in units of c; R is the 
dimensionless disk radius in units of the gravitational radius, 
Rg, integrated between characteristic inner and outer radii 
of -Rmncr and .Router; <f> ls the azimuthal angle of the disk to be 
integrated over (note that in the weak field approximation, 
a linear perturbation to the Special Relativity metric, the 
photons emitted at cf> make the same contribution as those 
emitted at n — <j>, and therefore the emission is integrated 
between — f and \, with the contribution to the emission 
from the back half of the disk accounted for by a factor of 2 
in the normalisation); i>o is the rest frequency and v is the 
observed frequency; and D = the Doppler factor, where 
v e is the emission frequency. 

In the weak field approximation, the Doppler factor is 

D= ^ Sd . (4) 

(1 + 7? 2 sin disk sin (j)) 
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Parameter 



Range 



Log? 



1 Flat continuum 

2 Broad H a shift with respect to narrow H a 

3 Broad H a central wavelength 

4 Broad H a Lorentzian width 

5 Broad H a intensity 

6 Narrow H a central wavelength 

7 Narrow H a Gaussian width 

8 Narrow H a intensity 

9 N II6550 /Narrow H a intensity ratio 

10 S II6718/S II6732 intensity ratio 

11 S II6732 /Narrow H a intensity ratio 

12 16302 /Narrow H a intensity ratio 

13 Second broad H a component shift with respect to nar- 
row H a 

14 Second broad H a component Gaussian width 

15 Second broad He* component intensity 

16 Disk intensity normalisation 

17 Disk shift with respect to 6564.61 A 

18 Sine of the disk angle 

19 Local velocity dispersion of the disk material 

20 Inner disk radius (-Ri nner ) 

21 Multiplication factor for disk outer radius 



-0.5 to 0.5 x 10- 9 Wm- 3 

0.985 to 1.015 Aref 

6540 to 6590 A 

20 to 1000 A 

0.05 to 20 x 10- 9 Wm" 3 

6540 to 6590 A 

1 to 30 A 

0.01 to 10 x 10~ 9 Wm- 3 

0.003 to 10 

0.2 to 3.0 

0.003 to 10 

0.003 to 10 

0.985 to 1.015 A ref 



20 to 1000 A 
0.05 to 20 x 10" 
10 to 10 5 
-100 to 100 A 
to 0.9994 
10" 3 to 10~ 2 c 
100 to 1000 R G 

2 to 100 Rinner 



1 Win- 



No 
No 
No 
Yes 
Yes 
No 
Yes 
Yes 
Yes 
No 
Yes 
Yes 
No 

Yes 
Yes 
Yes 
No 
No 
Yes 
Yes 
Yes 



Table 4. Prior ranges of input parameters for the emission models. Columns 1 and 2: Index and description of the 
free parameters; Column 3: Prior ranges of the parameter values; Column 4- Indication as to whether the coverage of 
the parameter space is logarithmic. Notes on units: lWm~ 3 = 10 _10 Wm~ 2 A~ 1 . A rc f is the wavelength of narrow Ho if 
present; the wavelength of the Lorentzian broad line component for models 5, 6 and 9; and the laboratory wavelength of 
H« in the case of model 12. 



Incorporating this into Equation ((3}, converting from fre- 
quency to wavelength, and simplifying: 



F v = K 



flout 



exp - 



XD 



-D 3 R" 



R — l 



(1 — sin 0diakcos 2 



dRd(j> 



(5) 



The parameters 7?i nn cr, Pouter, Omsk and Av were found 
from fitting to the observed line profile, while the normalisa- 
tion fixes eo-M 2 . The radial exponent of the disk emissivity 
q can also be fitted from the line shape; however, to reduce 
the number of free parameters and simplify the model, this 
parameter was fixed at a fiducial value of q — 3 for Ha 
l|Eracleous fc Halpernll20o3 ). This was a reasonable approx- 
imation to make, since the emission line flux re-radiated by 
the disk is proportional to the illuminating flux, which is 
predicted to vary as r~ 3 for a wide range of radii: the illu- 
minating flux falls as r~ 2 from the central source, and the 
flux falling per radius increment on the disk decreases asr -1 
due to geometric effects. 

The disk emission was calculated computationally us- 
ing a multi-dimensional Monte Carlo integration routine, 
gsLmontejuegas, from the GNU Scientific Library (GSL), 
for the grid of input parameters shown in Table [5] 

The wavelength coverage of the disk models is 6064 - 
7064 A, chosen as the Ha line emission was observed to 
be negligible outside this region for all quasars in the sub- 
sample. The emission models were calculated at a resolution 
of 10 A, which is adequate, as the resolution of the spectra 
themselves is 12 - 16 Apixel -1 , and the accretion disk emis- 
sion is smooth. The sine of the angle of the disk axis to the 



Parameter 


Range 


No. of 
points 


Log? 


Wavelength 


6064 to 7064 A 


101 


No 


Sine of the disk angle 


to 0.9994 


21 


No 


Velocity dispersion of disk 


10" 3 to 10" 2 c 


21 


Yes 


material 








Inner disk radius (-Rjnner) 


100 to 1000 R G 


21 


Yes 


Multiplication factor for 


2 to 100 -Rinncr 


26 


Yes 


disk outer radius 









Table 5. Ranges of input parameters for the accretion disk model. 
Column 1: Parameter description; Column 2: Range of parameter 
space covered by the models. The outer radius of the accretion 
disk is defined in terms of a multiplication factor for the inner 
radius; Column 3: Number of points in the array covering the 
range; Column 4- Indication as to whether the coverage of the 
parameter space is logarithmic. 



line of sight was allowed to vary across all of the possible 
parameter range, stopping just short of a value of unity (an 
exactly edge-on disk), since the disk has zero thickness in 
the model. The local velocity dispersion of the material in 
the disk covers the range 10 -3 to 10~ 2 c, which are typi- 
cal values for velocity dispersion in the BLR. The inner and 
outer disk radii were defined in terms of gravitational radii: 
the inner disk radius has a logarithmic range, and the outer 
disk radius was defined in terms of a multiplication factor 
for the inner radius. 

The disk emission models were checked both by eye 
and specially written automated routines, and those models 
found to have artifacts due to poor integration were recalcu- 



© 2008 RAS, MNRAS 000, [TJM] 



12 E. J. Down, S. Rawlings, D. S. Sivia and J. C. Baker 



lated using a larger number of steps in the integration rou- 
tine. During the Bayesian fitting process, the disk parame- 
ters (sine of the disk angle, sin(?disk, local velocity dispersion 
of the disk material, Av, and the inner and outer radii, Rumei 
and -Routor) were allowed to take continuous values over the 
prior ranges. The model for each set of parameters was cal- 
culated by interpolating between the sixteen disk emission 
templates which bracketed the required parameter values. 
This was a reasonable approximation, since the disk emis- 
sion varies slowly over each of the parameter ranges. 

It should be emphasised that the analysis of this paper 
is strongly dependent on the simplicity of the accretion disk 
model used, and the choice to fix the radial exponent of the 
disk emissivity to 3. A different disk morphology, such as 
an elliptical or warped disk, would alter the emission line 
profiles, and may affect the results. 



5 PARAMETER FITTING AND MODEL 
SELECTION 

5.1 Bayesian fitting 

The emission models were fit to the spectra following the 
Bayesian method, which uses a calculation of the likelihood 
of the recorded data arising, given a certain model, in or- 
der to find the probability distribution for each parameter. 
The Bayesian method is not discussed in detail here (see 
ISivia fc Skilungl l|2006l )). but it should be noted specifically 
that in the Bayesian context, the term "model" includes 
both the equation which describes the fit in terms of the 
free parameters, and the prior distributions of those param- 
eters. 

The Bayesian fitting was carried out using a "least 
squares" method. This folds in two important assumptions: 
first, that the prior is a fixed value over the entire range; 
second, that the noise on the data is well-approximated by 
a Gaussian distribution. There was limited prior knowledge 
available as to the values of the model parameters, so it was 
sensible to assign uniform priors over suitable ranges. The 
assumption that the noise on the data was Gaussian was a 
reasonable one; the noise on each data point was assumed 
to be Poisson-like sky emission noise, and since the data 
values were much larger than the error bars, this could be 
approximated as Gaussian noise. 

Sky emission dominates the errors in infrared spec- 
troscopy, so the error, Oi, on each data point was approxi- 
mated as the Poissonian noise resulting from sky emission, 
with an unknown scaling factor: 



Gi = QvSl, 



(6) 



where Q is the scaling factor, and S; are the values of the 
sky spectrum at each data point. The absolute values of 
the errors were not important to the analysis, and so the 
scaling factor Q was m arginalised. Following the method of 
ISivia fc Skilling (|2006T h the normalisation of the error bars 
was integrated out of the expression for the likelihood func- 
tion using a Jeffreys' prior on Q, which is uniform in log- 
arithmic space to encode ignorance as to the magnitude of 
the errors. 



5.2 BayeSys3 

T he Bayesian op timisation was carried out using BayeSys3 
bv lSkilling (20040, which explores the parameter space us- 
ing a range of Monte Carlo engines. BayeSys3 was called 
through a wrapper written by D. Sivia. For each model, a C 
program was used to calculate the likelihood function from 
each set of parameters provided by BayeSys3. 

BayeSys3 was initialised with an ensemble (i.e. the 
number of parallel explorations of the parameter space) of 
20; this was increased from an ensemble of 10 following the 
discovery of unstable results with smaller values (see Ap- 
pendix [X] for an overview of the tests performed to deter- 
mine the stability of the Bayesian fitting process) . All avail- 
able Monte Carlo exploration engines were switched on, to 
minimise the risk of probability density accumulating in lo- 
cal minima. The annealing rate, which controls the speed 
at which the simulation switches from exploring the entire 
parameter space to exploring the posterior parameter space, 
was set at what is suggested to be a rea sonably slow v alue 
for BayeSysS (0.1 in arbitrary units, see ISkillingl (12004 )). 



5.3 Model selection 

The Bayesian evidence is the probability of obtaining a cer- 
tain data set given a model, naturally weighted against mod- 
els with a larger prior parameter space. Models with unwar- 
ranted complexity are penalised in comparison to simpler 
models which fit the data equally well, according to Occam's 
Razor. The Bayesian evidence values from the twelve mod- 
els were compared for each quasar, to find the most likely 
model. Since the spectra have different resolutions, the evi- 
dence values of fits to different quasars are not compared. 
The Bayes factor Bab is written 



B A 



Prob(DjM A ) 



(7) 



Prob(D|M B )' 

where D is the data and M the model (|Trotta|[20oi ). The 
Bayes factor gives a statistical measure of the degree to 
which Model A has gained or lost support compared to 
Model B, given the d ata. The "Jeffreys' scale" shown in 
Table M (|Jeffrevslll939h provides an empirical scale for trans- 
lating the relative Bayesian evidence of two models into the 
more intuitive scale of odds, binning this into four bands of 
evidence: strong, moderate, weak and inconclusive. 

Based on the Jeffreys' scheme of Table [S] a quasar is 
considered to have strong evidence for the presence of a disk 
(denoted by "SD") if the model with the highest evidence 
is a disk model, and if no models without disks fall within 
one "Jeffreys' criterion" (A ln(Evidence) = 5) of the pre- 
ferred model. Moderate evidence for a disk (MD) and weak 
evidence for a disk (WD) are defined analogously, with the 
corresponding odds ratios. The quasar possesses a possible 
disk (PD) if the evidence is inconclusive, or if there is only 
weak or moderate evidence against the presence of a disk. 
The category of non-disk (ND) is assigned in cases where 
there is strong evidence against the presence of a disk, i.e. 
an accretion disk is apparently excluded by Jeffreys' crite- 



4 The BayeSys3 program and user guide are available at: 
htt p : / / www . inference . phy. cam .ac.uk /bayesy s / 
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|ln B AB | 


Odds radio 


Probability 


Strength of evidence 


<1.0 


<3:1 


< 0.750 


Inconclusive 


1.0 


~3:1 


0.750 


Weak evidence 


2.5 


~12:1 


0.923 


Moderate evidence 


5.0 


~150:1 


0.993 


Strong evidence 



Table 6. "Jeffreys' sc ale" for compa ring the strength of the evi- 
dence for two models (|jeffrevslll939h . Column 1 : Absolute value 
of the natural logarithm of the Bayes factor, which is a statistical 
measure of how Model A has increased or decreased in likelihood 
compared to Model B, given some data; Column 2: Odds ratio, 
an expression of the probability that Model A is true over Model 
B; Column 3: Probability of Model A compared to Model B; 
Column 4 '■ Empirical statement of the strength of the evidence. 



n(lnformation) 



Figure [3] illustrates the selection procedure with a plot 
of the natural logarithmic evidence versus the natural log- 
arithmic information for one quasar, although it should be 
noted that only the Bayesian evidence was used in the se- 
lection process. The quality of the data is a large factor in 
the model selection. Broadly speaking, high signal-to-noise 
spectra have very high ranges of evidence, as the difference 
between the best and worst fits is more apparent than for 
the lower signal-to-noise spectra, so the model selection pro- 
cedure is more conclusive. 

The Bayesian information is a measure of the ratio of 
the volume of the prior parameter space to the volume of 
the posterior parameter space. It therefore provides an in- 
dication of how much the data has increased knowledge of 
the parameter values, for a given model. There is a corre- 
lation between logarithmic evidence and logarithmic infor- 
mation for the high signal-to-noise cases: higher evidence is 
linked to better fits, which constrain the posterior parame- 
ters more tightly. For low signal-to-noise spectra, however, 
such as those of MRC0346-279 or MRC1217-209, there is no 
correlation apparent between the evidence and information, 
and one Jeffreys' criterion can encompass most of the mod- 
els, so it is impossible to discriminate reliably between them; 
in these cases, accretion disk emission may be present, but 
there is no evidence for it. 



5.4 Notes on individual quasars 

MRC0222-224: The model including the Lorentzian line, 
accretion disk and all narrow lines (Model 1) was selected 
with higher evidence by A ln(Evidence) = 6.4 than the next 
best model (Model 10), which includes the accretion disk 
and narrow lines only. 

MRC0327-241: The best-fit model includes the Lorentzian 
line, the accretion disk emission and narrow Ha (Model 3). 
Four other models are within the A ln(Evidence) = 5 bound 
of the selected model; however, all these models include disk 
emission, so this quasar has strong evidence for disk emis- 
sion. 

MRC0346-279: This spectrum suffers from low signal-to- 
noise, and so there is little evidence to discriminate between 
the models: all but one of the models fall within the Jef- 
freys' criterion of the best-fit model. Model 9, the fit with 
the broad Lorentzian and the broad Gaussian, was selected 
by the Bayesian evidence, which is the next simplest model 



Figure 3. Natural logarithmic evidence plotted against natu- 
ral logarithmic information for all model fits to the spectrum of 
MRC0450-221. The dotted line shows the evidence of the pre- 
ferred model, while the dashed line is plotted at an evidence dif- 
ference of one Jeffreys' criterion (i.e. A ln(Evidence) = 5) from 
the best fit. In this case, Model 1 is selected, with no other mod- 
els falling within the Jeffreys' criterion. 



after the single Lorentzian. The selected model has greater 
evidence by A ln(Evidence) = 0.81 than the next best fit. 
The low signal-to-noise ratio is likely to be the reason for 
the lack of fitted narrow emission lines in this spectrum. 
MRC0413-210: Model 1, which includes the Lorentzian 
broad line, the disk emission and all narrow lines, is pre- 
ferred, with A ln(Evidence) = 6.3 over the next best-fit 
model. 

MRC0413-296: This quasar is a special case, since there 
are obviously narrow lines present in the spectrum, but these 
did not constrain properly in the fit as the narrow Ho pa- 
rameters attempted to fit a broad component of emission. 
There appears to be more than two components of broad 
emission in this spectrum. From the evidence, the best model 
is Model 7, which has the broad Lorentzian, the broad Gaus- 
sian and all the narrow lines, although the narrow H a fitted 
a third broad component. A ln(Evidence) = 10.8 between 
the best fit, and the next best-fit model which includes a 
disk, so there is strong evidence against the presence of a 
disk. 

MRC0430-278: This spectrum has a reasonably low signal- 
to-noise ratio. The preferred model is Model 7, with the 
Lorentzian broad line plus the Gaussian broad line and all 
the narrow lines, though many models fall within the Jef- 
freys' criterion. The margin in logarithmic evidence of Model 
7 over Model 1, which includes the accretion disk emission, 
is only A ln(Evidence) = 2.2, so this is a "possible disk" 
quasar. 

MRC0437-244: Model 1, with the Lorentzian broad line, 
accretion disk emission and narrow lines, is preferred over 
the next best-fit model by A ln(Evidence) = 1.4, though the 
second best-fit model also includes accretion disk emission. 
There is strong evidence (A ln(Evidence) = 92.0) for the 
presence of a disk. 

MRC0450-221: The evidence is strong (A ln(Evidence) = 
40.1) for the presence of an accretion disk, with the selected 
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best-fit model being Model 1, including the accretion disk 
plus Lorentzian line and the full complement of narrow lines. 
MRC0549-213: Most of the models fall within the Jeffreys' 
criterion of the best-fit model, including all of the models 
with an accretion disk plus a standard BLR, so this quasar 
is a "possible disk" source. The preferred model for this low 
signal-to- noise ratio spectrum is Model 6, with one single 
Lorentzian line only. This is preferred over a disk emission 
model by A In (Evidence) = 2.3. 

MRC1019-227: This spectrum contains artifacts from 
poor sky subtraction which were excluded from the Bayesian 
fitting process, including the region surrounding narrow H a. 
As a result, the sub-set of models with the narrow Ha 
line only were not constrained, although the models with 
all narrow lines present could converge, since the narrow 
Ha parameters were fixed with reference to the other nar- 
row lines. The best-fit model is Model 1, with the broad 
Lorentzian, the accretion disk emission and all the narrow 
lines, though this is only favoured over Model 7 (which in- 
cludes the Lorentzian, the Gaussian and all the narrow lines) 
by A ln(Evidence) = 2.5, so this source is a "weak disk" can- 
didate. 

MRC1114-220: The selected model is Model 1, with the 
Lorentzian line plus accretion disk emission and the full com- 
plement of narrow lines. There is strong evidence for a disk, 
since the best-fit model without an accretion disk has lower 
evidence by A ln(Evidence) = 7.9. 

MRC1208-277: There is strong evidence supporting Model 
1, with all narrow lines and the Lorentzian line plus the 
accretion disk component. This has higher evidence than 
the best model without a disk by A ln(Evidence) = 34.8. 
MRC1217-209: Model 5, which includes the Lorentzian 
line and the accretion disk, is preferred, but the evidence 
for this model is weak. Most of the models fall within the 
Jeffreys' criterion of the best-fit model, with an evidence 
difference between Model 5 and the best-fit model without 
a disk of A ln(Evidence) = 1.4. The lack of fitted narrow 
lines is likely to be due to the poor signal-to-noise ratio of 
the spectrum. 

MRC1222-293: Model 1, which includes the Lorentzian 
line plus the accretion disk emission and the full set of nar- 
row lines, was selected as the best model with strong evi- 
dence ( A ln(Evidence) = 10.8). 

MRC1301-251: Model 1 (the Lorentzian and the disk emis- 
sion plus narrow lines), Model 7 (the Lorentzian and Gaus- 
sian plus narrow lines), and Model 10 (accretion disk emis- 
sion plus narrow lines), all have probabilities within the 
Jeffreys' criterion of each other and are clearly preferred 
above the other models. Model 7 is preferred by a margin of 
A ln(Evidence) = 0.5 over Model 1, so the results are incon- 
clusive. The wavelength shifts of the two broad components 
relative to the narrow H a line for Model 7 are —4200 km s _1 
for the Lorentzian line and +4500 km s _1 for the Gaussian 
line; these are within a plausible range for opposing outflows. 
It is clear, however, that more complex broad emission than 
a single broad line is required to fit this spectrum. 
MRC1349-265:The preferred fits are those with two broad 
components plus the full set of narrow lines. Of these, Model 
7 with the Lorentzian line plus Gaussian line and narrow 
lines is preferred, but only with weak to moderate evidence 
( A ln(Evidence) = 2.5) over Model 1, which includes emis- 



sion from an accretion disk in addition to the Lorentzian 
line and narrow lines. 

MRC1355-215: The preferred models for this spectrum are 
overwhelmingly those with two components of broad emis- 
sion and with all the narrow lines. Of these, Model 7, which 
includes the Lorentzian line and the Gaussian line in addi- 
tion to the narrow lines, is weakly preferred to Model 1, with 
the Lorentzian line and the accretion disk emission as well 
as the narrow lines, by A ln(Evidence) = 1.61. 
MRC1355-236: There is strong evidence 
(A In (Evidence) = 6.9) that Model 1 with the Lorentzian 
line, the disk emission, and all the narrow lines is preferred 
over the next best model. 

MRC1359-281: There is strong evidence 
(Aln(Evidence) = 6.6) for Model 1 with the Lorentzian 
line, accretion disk emission and all narrow lines. 



5.5 Correlations between the parameters 

Figure 2] shows the correlations between the posterior prob- 
ability for the accretion disk angle and the other disk fitting 
parameters of the selected model of the quasar MRC0450- 
221 (note that there are probability distributions from fits 
with five different random seeds marked on this plot). Simi- 
lar plots were examined for all fit parameters. There are no 
strong correlations of the local velocity dispersion of the disk 
material or the outer disk radius with any of the other disk 
parameters. There is a correlation such that when the sine of 
the disk axis angle increases, the inner radius increases and 
the wavelength shift of the disk component with respect to 
6564.61 A increases. 

In the model o f lChen fc Halpernl l| 19891 ) for an optically- 
emitting accretion disk, the emission line widths increase 
with decreasing inner disk radius. The inner disk radius, 
which is defined as a dimensionless quantity in units of 
the gravitational radius Rg, anticorrelates with the black 
hole mass. In the model for double-peaked lines, the line 
width therefore correlates positively with the black hole 
mass. Observationally, the widths of broad, low-ionisation 
emission line s are known to correlate with b lack hole mass 
ijVestergaardl l|2002h . lMc~Lure fc Jarvisl (|2002h l. which is pos- 
tulated to be due to virialisation of the emitting material. 

An increase in the disk angle causes the two peaks of the 
disk emission to move further apart, as the rotating material 
in the disk has a higher velocity component along the line 
of sight, as explained in Se ction | 4.5I This eff e ct ha s been 
observed observationally by iJarvis fc McLurd l|2006h . who 
found that using radio spectral indices as a proxy for source 
orientation, the sources at greater angles to the line of sight 
have larger broad- line widths. 

The correlation between the disk angle and the inner 
radius of the disk can therefore be interpreted in terms of the 
observed width of the broad emission; these two parameters 
act upon the emission line width in opposing senses, and so 
the fitted model is a trade-off between the two. 

The correlation of the disk angle with the shift of the 
disk component is caused by the stronger (blueward) line 
peak aligning with the strongest component in the spectrum, 
whereas the weaker (redward) peak of the line does not im- 
pose such a strong constraint on the fitting; the wavelength 
shift required to fit the line profile depends on the separation 
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Figure 4. A plot of the correlations between the sine of the 
disk angle (sinS^isk) an d the other accretion disk parameters for 
MRC0450-221: the inner disk radius (-Rinncr), the outer disk ra- 
dius (-Router), the local velocity dispersion of the disk material 
(Av), and the shift of the disk emission with respect to 6564.61 
A (Ai). There are distributions for five different random seeds 
plotted here: for all seeds, the probability distributions converge 
on similar regions (but see discussion in Appendix lAl . 



between the two peaks, which is most strongly dependent on 
the angle of the accretion disk to the line of sight. 

The disk angle has by far the greatest effect on the pro- 
file of the emission, since this parameter is most strongly 
correlated with the line width and the separation between 
the red and blue peaks of the line. The disk angle therefore 
constrains most strongly of all the disk fitting parameters, 
and so under the assumption of this given model, the double- 
peaked emission provides a reasonably robust method for 
measuring orientation. It should be emphasised that the 
model and input parameters chosen here (see Section l4~5j) af- 
fect the range of disk emission profiles available during the 
fitting process, and hence may affect both the disk angles 
fitted and whether there is Bayesian evidence for a disk. 



6 RESULTS AND DISCUSSION 
6.1 Fits to the emission spectra 

The best-evidence Bayesian fits to the quasar spectra and 
the residuals from these fits are shown in Figure where 
the fit is marked as a solid black line on the spectrum. In 
cases where disk models were selected, the disk component 
is shown by a dashed line, and the posterior distribution of 
the sine of the disk angle is shown in an inset plot. Each 
plot is labelled with the model with which it was fitted, and 
whether there is strong evidence for an accretion disk ac- 
cording to the Jeffreys' criterion (SD), moderate evidence 
(MD), weak evidence (WD), whether the results are incon- 
clusive but do not rule out a disk (PD), or whether there 



is strong evidence against the presence of a disk (ND) (see 
Table [S} . Figure [S] shows the best-fit models from the sub- 
set which include a disk, for those quasars whose selected 
best-fit models do not include an accretion disk. 

Table [7] summarises the results of the model selection 
process. Of the nineteen quasars, ten have strong evidence 
for disks according to Jeffreys' criterion; two have weak ev- 
idence for a disk; six have possible disks, which means that 
either the results are inconclusive, or that there is weak or 
moderate evidence against a disk; and in only one case is 
there strong evidence against the presence of a disk. 

The exceptional case, MRC0413-296, shows strong 
Bayesian evidence that there is no accretion disk emission in 
the spectrum. However, the best-fit model (Model 7, which 
includes the Lorentzian and Gaussian broad lines, in addi- 
tion to both narrow Ha and the forbidden narrow lines) 
did not fit in the expected manner: the width of the nar- 
row Ha emission was unconstrained in the fit, and this line 
was fitting to a third broad component. This quasar appears 
to be anomalous within the sub-sample in that two broad 
components are not sufficient to describe the emission, and 
therefore it is entirely plausible that this source requires a 
model not tested here, such as a three-component model in- 
cluding emission from an accretion disk, a BLR which gives 
rise to single-peaked lines, and an outflow. 

It is very notable that for all but two sources (PD 
cases MRC0346-279 and MRC0549-213), the selected mod- 
els include a complex BLR of more than one component. 
The most basic model of a single emitting region is not 
adequate to describe the complex profiles of the majority 
of these lines. In most of the cases, the preferred mod- 
els were the ones with all the narrow lines and either the 
Lorentzian broad line plus the accretion disk (Model 1), or 
the Lorentzian line plus the broad Gaussian line (Model 
7). In only four fits were models with less than the full 
complement of narrow lines preferred, and of these, three 
were spectra with low signal-to-noise ratios (MRC0327-241, 
MRC0346-279 and MRC1217-209). MRC0549-213 does ap- 
pear to be well-fitted with a single broad Lorentzian line, 
with the possible presence of weak forbidden narrow lines, 
but very weak or absent narrow H a. 

None of the models with the accretion disk emission 
only were selected; in most cases, these models have ex- 
tremely low evidence. There is certainly a component of the 
BLR which gives rise to single-peaked lines present. 

The fitted posterior probability distributions for the 
sine of the disk angle, shown as insets in Figures [5] and [6] are 
on the whole reasonably well-constrained with slightly asym- 
metric distributions, though in some cases, these are cut off 
by the zero-angle prior boundary. MRC1019-227 has a dou- 
ble peak in the posterior probability distribution, though the 
two peaks are closely spaced. MRC1217-209 has a poorly- 
constrained distribution, due to the low signal-to-noise ra- 
tio of the spectrum. The MRC0430-278 spectrum has a low 
signal-to-noise ratio, and the fitted posterior probability dis- 
tribution for the sine of the disk angle is therefore less tightly 
constrained in one region, although it is strongly peaked at 
low inclinations. MRC0549-213 has a posterior probability 
distribution for the sine of the disk angle which extends over 
the entire range. This source appears to have weak or ab- 
sent disk emission, as it is well-fit by a Lorentzian line, but 
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Figure 5. The best Bayesian fits to each quasar spectrum, and the residuals from the fit. The flux density scale of the 
residuals is the same as for the spectrum in most cases; for a few sources it is compressed, but the units remain the same. 
Wavelengths are rest-frame. The best fit is plotted as a solid black line. In cases where a disk is included in the selected 
model, the disk component is shown with a dashed line, and the posterior distribution of the sine of the disk angle is 
shown in an inset. In the inset plots, the y-axis shows the probability normalised to unity, and the x-axis covers the range 
< sin 6 < 1. The shaded area of the inset plots shows the la error bounds on the posterior distribution of the sine of the 
disk angle. Each plot is labelled with the index of the model with which it was fit, and with a code according to whether 
there is evidence for an accretion disk: SD = strong evidence for a disk; MD = moderate evidence for a disk; WD = weak 
evidence for a disk; PD = a possible disk (i.e. the evidence is inconclusive, or there is only weak or moderate evidence 
against the presence of a disk); ND = strong evidence against the presence of a disk. 



is classed as a PD source, as there is no strong evidence 
against a disk. 

The disks have fitted rotation axis angles between 1° 
and 48° to the line of sight; this range is consistent with 
the definition of quasars as being objects viewed within the 
opening angle of a dusty torus, where the opening angle is 



dependent on source luminosity but is ge nerally supposed to 
be roughly 45° for radio-luminous AGN ijLawrenc eT ll99ll ). 

Those sources with disk rotation axes at small angles 
to the line of sight in general have less strong evidence for 
a disk. The reason for this is that when the disk axis is at a 
small angle to the line of sight, the disk emission is not dis- 
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tinctive and double-peaked, but instead, rather similar to a 
Gaussian profile (see Section[53}. In these cases, there is lit- 
tle to distinguish the Lorentzian plus accretion disk model 
for the broad emission from the Lorentzian plus Gaussian 
model, except that the fit with the Gaussian has fewer pa- 
rameters, and is therefore likely to have a smaller prior pa- 
rameter space and hence be favoured by Occam's Razor. 
It is only those fits to sources with greater disk angles, or 
those with very high signal-to-noise ratio spectra, that make 
it possible to detect disk emission with high probability. 

Figure [7] shows the disk angles from all model fits for 
each quasar. In most cases, the angles are extremely stable. 
In cases where there are differences in the inferred angles, 
all angles from models within one Jeffreys' criterion of the 
best fit are consistent with each other. 



6.2 Relationships with projected radio source size 

The sine best-fit disk angles for the Molonglo quasars are 
plotted against the projected linear sizes of these sources 
in Fig ure [5] The source sizes are taken from iKapahi et al.l 
(1998), or in some cases, new sizes were obtained from ~ 
1.4 GHz MERLIN radio maps (Down et al., in prepara- 
tion) : this data is summarised in Table [5] The sample can 
be divided into three source types on the basis of previous 
studies: core-dominated sources, Compact Steep Spectrum 
(CSS) sources, and non-CSS FRIIs. Two of the 19 sources 
are core-dominated, which are sources whose jets are ori- 
ented at small angles to the line of sight, and therefore 
have core-to-lobe radio flux density ratios of greater than 
unity; examination of their radio-frequency spectral energy 
distributions (Down et al., in preparation) reveals that these 
quasars have been boosted into this sample by virtue of their 
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strong cores. The remainder of the sample is somewhat ar- 
bitrarily divided into six Compact Stee p Spectrum sources, 
with projected sizes of less than 25 kpc (|0 'Deaf l998) and 11 
non-CSS FRII sources, with projected sizes of greater than 
25 kpc. 

There is a highly tentative correlation with a Kendall 
r coefficient 0.44, probability of 65%, and signific ance of < 
la, calculated from ASURV (jLavallev et alJll992h . that for 
the eleven non-CSS FRII sources, the quasars with largest 
projected size have greater disk angles to the line of sight, 
as expected statistically from foreshortening if the accretion 
disk is perpendicular to the radio jet. The probability of a 
correlation is strengthened slightly to a Kendall r coefficient 
of 0.40 with probability of 74% and la significance if the six 
CSS sources are included in the sample. The weakness of the 



correlation is likely to be due to a large intrinsic scatter in 
the source size. 

There are four CSS sources with fitted disk angles less 
than 10° that appear to fall on the same relation as the non- 
CSS sources. It is probable that most or all of these are from 
the same population as the non-CSS sources, the only differ- 
ence being projected size. The remaining two CSS sources, 
MRC0222-224 and MRC1114-220, have strong evidence for 
disks inclined at angles greater than 30° to the line of sight. 

The radio map for MRC0222-224 (see iKapahi et ail 
l|l998l ) ) shows two radio lobes with no hint of a core (a pos- 
sible weak core is present between these lobes in a higher 
resolution MERLIN map, Down et al., in preparation); this 
quasar is consistent with an intrinsically small source viewed 
at a large angle to the accretion disk axis. The Balmer decre- 
ment for this quasar is estimated as Ha/ H/3 ~ 23, more 
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than a factor of two higher than for any of the other sources 
in this Molonglo sub-sample (Janssens et al., in prepara- 
tion). A simple interpretation is that this quasar is a young 
source, possibly surrounde d by a cocoon of dust which red- 
dens the optical emission (Ba ker et alJl20()3 ). 

The radi o map of MRC1114-220 l|de Silva et all 
lin preparation! ) reveals a strong radio core and single-sided 
jet, indicating that this source probably lies at a small angle 
to the line of sight, so that it is unlikely to be as intrinsi- 
cally small as it appears. Possible explanations are that the 
accretion disk and radio jet are misaligned in this source fol- 
lowing a merger event, or that the jet is precessing; however, 
this quasar merits further investigation. 

The two quasars excluded from the sub-sample, 
MRC0418-288 and MRC1256-243, have small projected 
sizes. MRC0418-288 is a CSS source with projected size 
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D < 16.6 kpc, so is likely to be an intrinsically small, red- 
dened source similar to MRC0222-224. MRC1256-243 has a 
larger size of D ~ 65 kpc and a high core-to-lobe flux den- 
sity ratio at 10 GHz of -Rioghz ~ 3, so is a core-dominated 
source: this is likely to have a small measured disk angle. 



6.3 Deprojected source sizes 

The source sizes were deprojected by dividing the appar- 
ent source sizes by the sine of the fitted disk angle, to 
compensate for simple geometric projection. This does not 
take into account the expansion of the source, but since the 
hotspots of jets are on ly expected to advance at ~ 0.1 c 
ijLongair fc Rilevlll979t ). this effect is small and is not con- 
sidered. The deprojected sizes are given in Table [8] and are 
plotted in Figure [3] 
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The distribution of deprojected source sizes found us- 
ing the fitted disk angles is a single-peaked distribution, with 
the main concentration of sources in the range 100 kpc - 1 
Mpc. Four of the CSS sources, MRC0222-224, MRC1114- 
220, MRC1349-265 and MRC1359-281, have deprojected 
sizes of less than 100 kpc. As the projected source sizes 
are defined as the distance between the cent res of the fur- 
thest separated radio components (following Ide Silva et al.l 
l|in preparation! ')'), the size may be underestimated if the 
whole source is not observed, i.e. in the case of MRC1114- 
220, only the core and approaching jet are visible, and so 
the size may be underestimated by a factor of ~ 2. There 
are hints that deprojection slightly tightens the distribution 
in source size as expected, but there is still a large scatter. 

The cumulative distribution of the deprojected source 
sizes found using the best-fit disk angles is shown in Figure 



1101 with expected distributions of linear sizes, assuming that 
the hotspots of the lobes are expanding at a constant rate, 
shown for comparison. If the radio jet and accretion disk 
axes coincide, the calculated deprojected sizes are broadly 
consistent with a constant expansion of the heads of the 
jets up to ~ 1 Mpc, although there seems to be an excess 
of small sources. The distribution drops off at sizes greater 
than ~ 1 Mpc, which can be explained by the duty cycle 
of the quasars: the number of sources larger than a cer- 
tain cut-off value will be depleted as they become quiescent 
ijBird et al . 2008). The distribution tails off more gradually 
than the highly simplified model, which can be explained 
by variation in hotspot advance speeds between quasars, 
some sources having an unusually long duty cycle, or by the 
largest sources being in especially low-density environments. 
The excluded sources, MRC0418-288 and MRC1256-243, are 
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Figure 6. The best Bayesian fits which include an accretion disk component, for each quasar for which the highest evidence 
fit did not include a disk. Residuals are plotted below each fit. The flux density scale of the residuals is the same as for the 
spectrum in most cases; for a few sources it is compressed, but the units remain the same. Wavelengths are rest-frame. The 
best fit is plotted as a solid black line. The disk contribution is shown with a dashed line, and the posterior distribution of 
the sine of the disk angle is shown in an inset. In the inset plots, the y-axis shows the probability normalised to unity, and 
the x-axis covers the range < sin 6* < 1. The shaded area of the inset plots shows the la error bounds on the posterior 
distribution of the sine of the disk angle. Each plot is labelled with the index of the model with which it was fit, and with 
a code according to whether there is evidence for a disk: SD = strong evidence for a disk; MD = moderate evidence for 
a disk; WD = weak evidence for a disk; PD = a possible disk (i.e. the evidence is inconclusive, or there is only weak or 
moderate evidence against the presence of a disk); ND = strong evidence against the presence of a disk. 



predicted to have very small and very large deprojected sizes 
respectively, so will not affect the overall distribution much. 



6.4 Relationships with radio luminosity 

The sine of the best-fit disk angles are plotted against the 
178 MHz radio luminosity in Figure 1111 There is a correla- 
tion between these quantities, with a Kendall r coefficient of 
0.60, probability of 93% and 1.8cr significance. The measure- 
ment of the disk angle is entirely independent of the radio 
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luminosity, so the fact that these parameters are correlated 
provid es direct, albeit weak, evidence for the receding torus 
model |Lawrencell 199 ll ) , if the accretion disk axes align with 
the radio jets. According to this theory, sources with high ra- 
dio luminosit y have larger t orus opening angles due to dust 
sublimation (jSimpsonll 19981 ). and therefore higher luminos- 
ity sources appear as quasars, rather than radio galaxies, up 
to a greater viewing angle (the critical angle). The predic- 
tion from this model is therefore that more luminous quasars 
have jet angles ranging up to a higher cut-off value, and if 
the disk angle is assumed to match the jet angle, the calcu- 
lated disk angles support this prediction. 

A ca lculation of the toru s opening angle, following the 
model of IWillott et ail (2000), normalised by a critical an- 
gle of 45° at logLi78MHz= 27, and modified by a minimum 
quasar fraction of 10% which dominates the opening angle 



at low radio luminosities (e.g. Vardoulaki et all (120081 )1. is 
also marked on Figure ITTI The measured disk angles all fall 
within the calculated envelope of opening angles, and so as- 
suming that the obscuring tori are aligned with the accretion 
disks, these results are in accordance with the receding torus 
model. 



The missing source MRC1256-243 has luminosity 
log Li78MHz ~ 27.4, which is consistent with a low disk an- 
gle. MRC0418-288 has luminosity log Li 78 mhz ~ 26.8, which 
is a relatively low radio luminosity for the predicted large 
disk angle; however, an angle of ~ 30° would still be consis- 
tent with the receding torus scheme. 
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Table 7. Summary of the model fitting results. Columns 1 and 2: Quasar index and MRC name; Column 3: Index of 
the best-fit model. If the selected model does not include an accretion disk, the index of the highest-evidence model which 
includes an accretion disk is shown in brackets. Column 4- Indication of whether there is evidence for a disk according to 
the odds ratio, with codes as follows: SD = strong evidence for a disk; MD = moderate evidence for a disk; WD = weak 
evidence for a disk; PD = a possible disk, i.e. the best fit is a non-disk model, but there is only inconclusive, weak or 
moderate evidence for this; ND = strong evidence against the presence of a disk; Columns 5 and 6: Natural logarithmic 
evidence difference between the best-fit model, in cases where this includes an accretion disk, and the second best model 
{Column 5. A In E^,), and in cases where the second best fit also includes a disk, between the best- fit model and the best-fit 
model without a disk (Column 6, AlnE" d ); Columns 7 and 8: Natural logarithmic evidence difference between the best-fit 
model, in cases where this does not include an accretion disk, and the second best model (Column 7, Aln.E 2 n ), and in 
cases where the second best fit also does not include a disk, between the best-fit model and the best-fit model with a disk 
(Column S, Aln Ef n ); Column 9: Best-fit angle of the disk and the error in degrees. 



6.5 Distribution of angles 

The expected distribution of jet angles (9) can be mod- 
elled for the Molonglo quasar sample, including a Doppler- 
boosted core component. From Bayes' Theorem 

Prob(6»|L tot > L min ) oc Prob(L to t > L min |6») x Prob((9) (8) 

where L to t is the total luminosity of the quasar and L m i n is 
the limiting luminosity of the survey, and 



Prob(<9) = 27rsin6>d<9 



0) 



The proportionality is required because there is no solid in- 
formation about the overall number of sources for which 

Ltot > ^min- 

Marginalising the lobe luminosity, then 



Prob(L to t > -kminl^) oc / Prob(L t ot > L min \9, L lobc ) 

Jo 

xProb(Liobc)dLi ob c • (10) 

From IWillott et al.l (|200ll ). the radio luminosity func- 
tion of high luminosity AGN is 



p{L) = Pn 



L 



(11) 



where p is the comoving space density of sources in logarith- 
mic luminosity space and a ~ 2.3. Since this relation was 
found for luminosities of 151 MHz and 178 MHz, then the 



total luminosity is low enough in frequency to be approxi- 
mated as the lobe luminosity. Then 



Prob (Liobc) oc L lob " e x 



oc 



dlog Liobc 

Llobc 



(12) 



Defining the core-to-lo be flux ratio R in the same way as 
IJackson fc Wail (|l999h . then 



R — R c —^ 



7 2 V(l + I3cos9) 2 (l-/3cos(9) 2 



(13) 



where R c is some fid ucial value of R, and is found by 
IJackson fc Walll (| 19991 ) to be R c ~ 0.01 for FRII sources; 
7 is the Lorentz factor and (3 is the velocity of the jet in 
units of c. The total flux density, Stot, is then defined as 

Stot = S*iobc(l + i?) , (14) 

and since flux density is proportional to luminosity, then 

Ltot = Li obe (l + R) . (15) 

Now for each value of Liobe and 9, the total luminosity 
is uniquely defined, and Prob(L to t > L m i n |0, Liobc) becomes 
simply or 1. The outcome of this is that for a given 9, there 
is one limiting lobe flux Ln m (6 l ) above which the probability 
of detecting a source is unity, and below which it is zero, and 
this simply changes the limits on the integration so that, 
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Figure 7. Disk angles for model fits with an accretion disk. Circu- 
lar symbols indicate evidence for a disk, square symbols indicate 
a possible disk (inconclusive evidence, or weak to moderate ev- 
idence against a disk), while triangular symbols indicate strong 
evidence against the presence of a disk. White points mark fits 
which fall within one Jeffreys' criterion of the best fit, and these 
can be seen in all cases to be consistent with the best-fit disk 
angle. The remainder of the fitted disk angles are plotted in pale 
grey. 



substituting equations ©, JTUJ and (|12|) into equation (JSJ 
we find 



Prob(6<|Ltot > £min) oc 

This integrates to 

Prob(6»jLtot > Lmin) 



poo 
J L lim (8 



d-Liobe x sinOdfl . (16) 



i — a I oc 



OC Llim(0)" 



m(") X Slll( 

x sin6> . 



(17) 



Figure [121 shows the cumulative probability distribution 
for the disk angles, together with the theoretical distribu- 
tions of jet angle for different values of the Lorentz factor 7. 
If the accretion disk axis and radio jets align (i.e. if the disk 
is perpendicular to the jet), then the distribution of disk an- 
gles is most consistent with a value of the jet Lorentz factor 
around 7 ~ 20. It should be noted that MRC0418-288 was 
excluded from this sub-sample of quasars on the basis of op- 
tical faintness. This source has a high probability of being 
at a large angle to the line of sight, since it may well be a 
reddened source in which the dusty torus is bl ocking some 
sight lines to the optically bright nucleus (e.g. iBaker et al.l 
(2002)). MRC1256-243 is a core-dominated source, which is 
likely to have a very small disk angle. An extra quasar added 
to the each end of the disk angle distribution function would 
make it more in agreement with the modelled angular dis- 
tribution for 7 = 20. 

This Lorentz factor differs from the 7 ~ 3 measured 
from the 3C sample. This is explicable if 7 is dependent on 
angle, such that 7 ~ 20 close to the line of sight, when view- 
ing the quasars down the ax es of their jets, an d 7 ~ 3 when 
viewed at larger angles (e.g. lHardcastlel l|2006l )). The angu- 
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Figure 8. Sine best-fit disk angles versus the projected linear 
size of the quasar in kpc. The disk angle error bars include both 
the errors returned from the individual fits and a systematic error 
calculated from the variation in results from using different ran- 
dom seeds to initialise the Monte Carlo engines, added in quadra- 
ture, but do not take into account any error due to model se- 
lection. The non-CSS FRII sources are plotted with circles, and 
CSS sources are plotted with squares. The two sources plotted 
with crosses are core-dominated sources without measured pro- 
jected sizes. The black points have weak to strong evidence for 
a disk, grey points have possible disks, and the white point has 
Bayesian evidence against the presence of a disk. Index numbers 
are matched to source names in Table [T] The black line shows 
the best linear fit to all sources excluding the core-dominated 
sources. There is a highly tentative correlation (Kendall r coeffi- 
cient of 0.40, significant at the la level) between the sine of the 
disk angle and the projected linear size for all sources excluding 
the core-dominated sources, and for the 11 non-CSS FRII sources 
(Kendall r coefficient of 0.44). The dotted, dot-dot-dot-dashed, 
dot-dashed, dashed and long-dashed lines show loci of constant 
deprojected size, for 10 kpc, 25 kpc, 100 kpc, 1 Mpc and 5 Mpc 
respectively. 



lar distribution of the Molonglo quasars, which are selected 
at intermediate radio frequency, is dominated by Doppler 
boosting, whereas this is a small effect for the 3CRR sample. 
The distribution of the fitted Lorentz factor for the 3CRR 
sample is centred around 7 ~ 3, though with a slight asym- 
metry biased towards higher values of 7, whereas the Molon- 
glo sample is predicted to have more probability of higher 
7 factors. It is not sufficient to model the jet with a single 
Lorentz factor. The luminosity - jet angle relation is also ex- 
pected to have some scatter in 7 due to intrinsic differences 
in the sources. 

The distribution in fitted disk angle levels off at 48°. 
This is consistent with the opening angle which might be 
expected for a sample of powerful FRII sources. The fact 
that no quasars are observed with disk angles larger than 
48° adds to the evidence for the unification scheme which 
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Table 8. Summary of projected and deprojected source sizes. 
Columns 1 and 2: Quasar index and MRC name; Column 3: 
Source type - CD indicates sources with core-to-lobe radio flux 
ratio at fO GHz of greater than 1, CSS indicates sources with 
ctopt > 0.5 and projected source sizes < 25 kpc, FRII indicates 
non-CSS FRII, which encompasses all other sources; Column 4 '■ 
Projected source sizes in kpc; Column 5: Origin of projected 
source size measur ement - 1 = New me asurement (Down et al., in 
preparation), 2 = iKapahi et al.l ljl998T ): Column 6: Deprojected 
source sizes in kpc calculated using the fitted disk angles. 



suggests radio galaxies and quasars are the same objects 
viewed at decreasing angles to the line of sight. 

6.6 Velocity shifts 

6.6.1 Velocity shift measurements 

The prior ranges on the positions of the single-peaked broad 
line components allowed shifts of Az — ±0.015 compared 
to the fitted positions of the narrow lines, corresponding 
to line shifts of ~ 4500 km s _1 . If no narrow lines were 
fitted, the shift of the Gaussian line was measured relative 
to the Lorentzian line. The accretion disk was allowed to 
shift to a similar degree, AA = ±100 A with respect to the 
laboratory wavelength. The fitted velocity shifts from the 
best-fit models are given in Table [9] 

In seven cases, the measured velocity shifts cannot be 
trusted. This occurs in cases where the fitted parameters 
of either the broad Lorentzian line or narrow Ha are not 
constrained by the prior range. There are four potential rea- 
sons for this. (1) The velocity shift of the broad Lorentzian 
is unconstrained by the fit; this component is fitting to a 
broadband bump in the continuum left by imperfect con- 
tinuum subtraction, rather than to the emission line profile, 
so this fitted value is not correct. This affects the quasars 
MRC0327-224 and MRC1301-251. (2) The broad Lorentzian 
is clearly fitting to an additional narrow line component, and 
has unconstrained width as a result. This affects MRC 1208- 
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Figure 9. The solid line shows source sizes deprojected with the 
fitted disk angles (see Table [HJ. The projected sizes are plotted 
with a dotted line. Note that the two core-dominated sources have 
no measured size and are therefore excluded. 




Deprojected source size D / kpc 

Figure 10. Cumulative distribution of source sizes deprojected 
using the best-fit disk angles (solid line), with the projected source 
sizes shown for comparison (dotted line). The curves show the pre- 
dicted distribution if the heads of the sources are expanding at 
a constant rate; the dot-dashed line is normalised at the depro- 
jected size of the largest source in the sample, the dashed line is 
normalised at 1 Mpc. Note that the two core-dominated sources 
have no measured size, and are boosted into this sample by strong 
core emission, and are therefore excluded from this plot. 



277. (3) Narrow H a is fitting to an extra broad compo- 
nent, and therefore the fitted position of this line is not 
reliable. The relative shifts of the broad components to each 
other will be of the right magnitude so long as the extra 
broad component fitted by narrow H a is small. This affects 
MRC0413-296, MRC0430-278 and MRC1222-293. (4) In the 
case of MRC0222-224, there is some degeneracy in the fit of 
the broad Lorentzian, such that the shift of this line is not 
properly constrained. 
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Table 9. Velocity shifts of the broad line components. Columns 1 and 2: Quasar index and MRC name; Column 3: Best- 
fit disk model; Column 4-'- Disk emission shift relative to narrow He? in km s _1 for the best-fit disk model; Column 5: 
Lorentzian line shift relative to narrow Ho in km s — 1 for the best-fit disk model; Column 6: Best- fit model in cases where 
this does not include a disk; Column 7: Broad Lorentzian shift relative to narrow He* in km s - i, for the best-fit model in 
cases where this does not include a disk; Column 8: Broad Gaussian shift relative to narrow Ha in km 8 , for the best-fit 
model in cases where this does not include a disk. * No narrow H a was fitted in this model, and so the shift of the disk 
emission relative to the broad Lorentzian is given in km s — 1 . ' No narrow Ho was fitted in this model, and so the shift of 
the broad Gaussian relative to the broad Lorentzian is given in km s . 



6.6.2 Discussion of the velocity shifts 

Figure [T3] shows the shifts of the broad emission components 
for the best-fit models which include an accretion disk. Note 
that these models all include a Lorentzian line in addition 
to the disk emission. Three of the sources were best fitted 
with models which did not include narrow lines, and hence 
do not have shifts measured relative to the NLR; none of 
these quasars have published redshifts from high-resolution 
observations of narrow lines, so these were excluded. The dis- 
tribution of fitted disk shifts is relatively evenly distributed 
between -2000 km s _1 and +2000 km s -1 , with the excep- 
tion of MRC0437-244 with a disk shift of -2660 km s" 1 rel- 
ative to the narrow lines. This distribution does not change 
significantly if only the reliable shifts are considered. 

The fitted Lorentzian line shifts are distributed in three 
groups. The largest group have relatively modest velocity 
offsets to the NLR, with 9/16 shifts within ±500 km s" 1 
of the narrow line position, and a small tail out to -2000 
km s . There is one source only with a Lorentzian line 
shift relative to the narrow lines of more than +500 km s -1 : 
MRC0413-296 has a fitted shift of +2750 km s~\ but this is 
not reliable, since the narrow lines were fitting to an extra 
broad component in this spectrum and hence the Ha line 
centre used for reference was not reliably fitted. There are 
also three sources with blueshifts of -3500 - -4500 km s" 1 . Of 
these three quasars, two have Lorentzian lines which are fit- 
ting to a component of the continuum rather than the broad 
H a emission, and the other has degeneracy in the fitted val- 
ues of the Lorentzian line, so the shifts for these sources are 



not trustworthy. Therefore, all dependable Lorentzian line 
shifts fall in the central group, with moderate redshifts with 
respect to the NLR, or moderate to high blueshifts with 
respect to the NLR. 

Table [10] gives the average shifts of the broad line com- 
ponents and their standard errors. The spread of shifts is 
very wide, but when only the reliable velocity shifts are con- 
sidered, both the accretion disk and the Lorentzian line have 
average blueshifts of order -200 km s . Taken overall, there 
is a trend within this sample of powerful quasars for the 
broad emission to be blueshifted, in agreement with the ma- 
jority of the literature (see Section fl.3|l . It should be borne 
in mind that the outflows and inflows may be symmetric, 
but if dust in the plane of the accretion disk obscures the 
far-side emission from view, the shifts measured are for the 
near-side quasar emission. 

Separating the sample into 10 low-z (z < 1) and 6 high- 
z (z > 1.5) sources and considering only the trustworthy 
shifts (6 and 3 sources respectively for low-z and high-z), 
the average Lorentzian line shift is -380 km s~i for the low- 
z sources and +40 km s _1 for the high-z sources, consistent 
with this emission originating close to the systemic velocity. 
The average accretion disk shift for the trustworthy sources 
is -890 km s for the low-z sources and +620 km s for 
the high-z sources, so the disk emission is blueshifted with 
respect to systemic velocity at low-z and redshifted at high- 
z. 

Figure [141 shows the broad-line velocity shifts from the 
best-fit disk models plotted against quasar redshift. The red- 
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All shifts 






Trusted shifts only 




Source of emission 


Average shift 


Standard error 


No. of sources 


Average shift 


Standard error 


No. of sources 


(and z range) 


(km s L ) 


(km s x ) 




(km s 


(km s 1 ) 




Accretion disk (all) 


-400 


310 


16 


-180 


470 


9 


Accretion disk (low-z) 


-920 


320 


10 


-890 


460 


6 


Accretion disk (high-z) 


460 


440 


6 


620 


260 


3 


Lorentzian line (all) 


-880 


470 


16 


-240 


230 


9 


Lorentzian line (low-z) 


-1210 


520 


10 


-380 


320 


6 


Lorentzian shift (high-z) 


-330 


890 


6 


40 


150 


3 



Table 10. Summary of the mean averages and standard errors of the shifts of the broad emission components, for best-fit 
models which include an accretion disk. Shifts arc given for both the entire sample, and only those sources whose shifts are 
reliable; and are given for sub-samples divided by redshift: the sample as a whole, the low-z fraction (z < 1), and the high-z 
fraction (z > 1.5). The number of sources in each sub-sample is shown. 
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Figure 11. Sine best-fit disk angles versus logarithmic 178 MHz 
luminosity. The disk angle error bars include both the errors re- 
turned from the individual fits and a systematic error calculated 
from the variation in results from using different random seeds 
to initialise the Monte Carlo engines, added in quadrature, but 
do not take into account any error due to model selection. The 
non-CSS FRII sources and core-dominated sources are plotted as 
circles, and CSS sources are plotted as squares. The black points 
have weak to strong evidence for a disk, the grey points have 
possible disks (there is inconclusive evidence, or weak to mod- 
erate evidence against a disk) and the white points have strong 
evidence against a disk. Index numbers are matched to source 
names in Table [l] The best-fit line (minimising x 2 ) is shown in 
black. The dashed li ne shows the r a dio lu minosity-dependent crit- 
ical angle (following Iwillott et all fcOOCh ) with fiducial angle 45° 
at logLi78MHz= 27, and minimum quasar fraction 10%. 



shift and the shift of the accretion disk emission with re- 
spect to the narrow lines are correlated with a Kendall r 
coefficient of 0.850, significant at the 2<r level (98% proba- 
bility of a correlation). This correlation may arise from the 
dependence of both these quantities on optical luminosity 
(Wills 1980) or radio luminosity. The relation between AGN 
redshift and luminosity is well-known, and occurs because 
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Figure 12. The cumulative distribution of the best-fit disk angles 
(solid black line). The curves s how the theoretical dist ributions 
of jet angle, modelled from the Ja ckson fc Walll Jl999t) relation 
and normalised up to an angle of 45°, for different values of the 
Lorentz factor: the dashed line shows the expected cumulative 
distribution for 7 = 3, the dot-dashed line for 7 = 8.5, and the 
dotted line for 7 = 20. The vertical dot-dot-dot-dashed line indi- 
cates a disk angle of 45° . 



sources at greater distances are more powerful. This effect, 
the Malmquist Bias, is illustrated in Table flTl which shows 
that the sources with z > 1.5 are more luminous than the 
z < 1 sources at radio frequencies, and in terms of opti- 
cal continuum and emission lines. There is no correlation 
apparent between the quasar redshift and the shift of the 
fitted Lorentzian line with respect to the narrow lines, with 
a Kendall r of 0.133 (significance of < 0.5<r, 28% probabil- 
ity of a correlation). It appears that the shifts between the 
single-peaked broad emission line and the NLR are modest. 

Large-scale emission is most commonly seen alongside 
the radio jets or w i thin t he ionisation cone of the source, 
e.g. iNesvadba et all (|200ct ) . If it is the case that the material 
flows alongside the jet, the greatest velocity shifts should be 
observed for sources at small angles to the line of sight. No 
correlation is found between the magnitude of the velocity 
shift of the disk with respect to the NLR and the sine of the 
disk angle to the line of sight (Kendall r coefficient -0.05, 
significance of < O.lcr, probability of correlation 11%); the 
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Figure 14. Accretion disk shifts (circles) and Lorentzian line shifts (squares) with respect to narrow H a for the best- 
fit models which include a disk, plotted against the source redshift. Black points indicate sources in which the shifts are 
reliable, and grey points indicate sources in which the shifts are not trusted (see Section |6. 6, It . Dotted lines join the disk 
and Lorentzian shifts from the same source, where the error bars do not overlap. The lower limits of the error bars for 
sources 2, 15 and 1 which are not shown on the plot are -9700 km s —1 , -7600 km s — 1 and -9400 km s — 1 respectively. 





Z < 1 




z > 1.5 




All z 






Av. 


Err. 
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Err. 


Av. 


Err. 


-^lslMHz 


26.9 


0.1 


27.8 


0.1 


27.3 


0.1 




-24.8 


0.3 


-25.7 


0.4 


-25.1 


0.3 


^Nar H a 


0.135 


0.030 


0.541 


0.221 


0.285 


0.095 


-^Brd H a 


3.93 


1.05 


19.54 


6.49 
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Table 11. Average luminosities for the quasar sub-sample di- 
vided by redshift into a low-2 fraction (z < 1), and a high-z 
fraction (z > 1.5), and for the combined sample. Row 1: Log- 
arithmic 151 MHz luminosity; Row 2: Optical bj magnitudes 
from the UK Schmidt Ij aJ plates, where 6j = B - 0.23(5 - V) 
( Bahcall fc Soneirall 19801) ; Row 3: Luminosities of the narrow Ha 
lines in units of 10 36 W sr _1 ; Row 4'. Luminosities of the broad 
Ha lines in units of 10 36 W sr — 1 . 



velocity shifts are plotted against the disk angle in Figure 
1151 However, the sources at > 30° show velocity shifts of less 
than ±1000 km s _1 (note that these are not reliable shifts), 
with the exception of MRC1114-220, which is a young radio 
source, and hence is more likely to have a misaligned disk. 
The narrow line gas may in fact be entrained along the jet , 
but moving in a transverse direction: iNesvadba et alJ (|2006l ) 
and lNesvadba et all l|2008l ) suggest from integral field stud- 



ies of a small sample of high-redshift radio galaxies that 
the ionised gas is being accelerated away from the nucleus 
of the galaxy by expanding cocoons of hot gas around the 
jets. There are also indications that the distribution of gas 
surrounding CSS sources does not mirror that seen in the 
larger sources, with more dust and gas seen in edge-on CSS 
sources, and so it is probable that the velocity distribution 
of the material in these young sources will also differ. 

There is no evidence for a correlation between the mag- 
nitude of the disk shift with respect to the NLR and the 
width of the narrow lines (deconvolved from instrumental 
effects), with Kendall r of -0.283, significance < la, a prob- 
ability for anticorrelation of 56%. The narrow lines would be 
expected to have broader profiles for greater velocity shifts 
if the shift between the disk emission and the narrow line 
emission is due to the narrow lines forming in a moving mass 
of gas. 

The measured velocity shifts for the broad components 
relative to the narrow lines imply that the single-peaked 
broad lines arise from a region with a similar velocity to the 
NLR, while the disk is redshifted relative to these regions 
for higher-z, and hence brighter, sources, and blueshifted in 
relation to these for lower-z sources. The accretion disk is 
expected to be fixed at the systemic redshift, and therefore 
these results can be interpreted as the single-peaked broad 
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Figure 13. Distribution of line shifts for the best-fit models 
which include an accretion disk. The disk shifts with respect 
to narrow H a are shown as a grey-shaded histogram, and the 
Lorentzian line shifts with respect to narrow H a are shown as 
a solid-line histogram. The subsets of these with reliable mea- 
surements are marked with dark-grey shading and dot shading 
respectively. Note that three sources are excluded as no narrow 
lines were fit for these spectra. 
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Figure 15. Best-fit accretion disk angle plotted against the shift 
of the disk emission with respect to narrow H a for non-CSS 
sources (circles) and CSS sources (squares). The points with reli- 
able shifts are plotted in black; those sources whose shifts are not 
trusted are plotted in grey. 



lines and narrow lines arising from outflowing material in 
optically- luminous sources, and infalling material in less lu- 
minous sources. This could be the result of a quasar wind 
model of AGN feedback, in which powerful qua sars fuel out- 
flows which sweep the gas out of the galaxy (|Silk fc Reesl 



Il998t) : the outflows are expected to scale with the quasar 
luminosity. 



7 EFFECTS OF MODEL ASSUMPTIONS 

7.1 Models with the NLR and accretion disk 
fixed at the same redshift 

If the narrow lines are formed at the systemic redshift, as 
much of the literature indicates, then there should be mini- 
mal velocity shift between the accretion disk and the NLR. 
In order to test the effect of this alternative assumption, a 
second set of models was fitted, with the disk emission fixed 
at the same redshift as the narrow lines. Aside from the ex- 
clusion of one parameter describing the shift between these 
two emission components, the models are identical, and in 
the cases of the models with no narrow lines, the disk is 
allowed to shift as usual; this change therefore affects only 
four models (1, 3, 10 and 11). 



7.2 Notes on individual quasars 

MRC0222-224: Model 1 was selected with higher evidence 
by A ln(Evidence) = 6.0 than Model 10, the next best 
model. 

MRC0327-241: The best-fit model is Model 1. Models 3, 
5 and 10 are also within the A ln(Evidence) = 5 bound of 
this model; however, all these models include disk emission, 
so this quasar has strong evidence for a disk. 
MRC0346-279: This low signal-to-noise spectrum does 
not differentiate well between the models. Model 9 was se- 
lected by the Bayesian evidence, with higher evidence by 
A ln(Evidence) = 0.80 over the next best fit. 
MRC0413-210: Model 1 was selected, with 
A ln(Evidence) = 7.7 over the next best model. 
MRC0413-296: As before, although the spectrum of this 
quasar clearly contains narrow lines, the narrow line param- 
eters instead fitted broad components of emission. The best 
model is Model 7, with evidence over the next best-fit model 
of A ln(Evidence) = 12.1. This quasar has strong evidence 
against the presence of a disk. 

MRC0430-278: The preferred model for this spectrum is 
Model 7, though many models fall within the Jeffreys' cri- 
terion. The margin in logarithmic evidence of Model 7 over 
Model 3, the best-fit model which includes disk emission, 
is only A ln(Evidence) = 2.3, so this is a "possible disk" 
quasar. 

MRC0437-244: Model 5, including a disk, was selected 
with strong evidence. 

MRC0450-221: Model 1, including a disk, was selected 
with strong evidence. The fitted disk angle is not within the 
errors of the fitted angle from the first analysis: 20° ± \ from 
this analysis compared to 13° ± i for the case in which the 
accretion disk is allowed to shift. 

MRC0549-213: Many of the models fall within the Jef- 
freys' criterion of the best-fit model, including three of the 
disk models, so this quasar is a "possible disk" source. The 
selected model was Model 6, but this is preferred over the 
best-fit disk model, Model 1, by only A ln(Evidence) = 1.7. 
MRCT019-227: Model 1 is preferred by only 
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A ln(Evidence) = 1.5 over Model 7, so this source has 
a "weak disk" classification. 

MRC1114-220: Model 1 was selected by A In (Evidence) = 
6.3 over Model 7, the best-fit model without a disk, so there 
is strong evidence for a disk. The fitted angle, 15° ± 4, is 
not in agreement with the angle found from the analysis in 
which the disk is allowed to shift with respect to the narrow 
lines (35° ±\ 2 ). 

MRC1208-277: Model 7 was selected with evidence differ- 
ence A ln(Evidence) = 86.3 over Model 3, which contains a 
disk, so this has strong evidence against a disk. The fitted 
angle is 46° ± 3, which is strongly in disagreement with the 
case in which the disk is allowed to shift with respect to the 
NLR (15° ±\). 

MRC1217-209: Model 5 was selected, but many of the 
models fall within the Jeffreys' criterion of this model. There 
is an evidence difference between Model 5 and Model 9, the 
best-fit model without a disk, of A ln(Evidence) = 1.4. 
MRC1222-293: Model 7 was selected as the best-fit model, 
with evidence difference of A ln(Evidence) = 50.4 over 
Model 1, which is the highest-evidence model including an 
accretion disk. The fitted angle is 19° ± 3, which is in dis- 
agreement with the angle found from the analysis is which 
the disk is allowed a velocity shift with respect to the NLR 
(4°±f). 

MRC1301-251: Model 7 is preferred by a margin of 
A ln(Evidence) = 0.05 over Model 10, the highest-evidence 
model including a disk, and by a margin of A ln(Evidence) = 
0.8 over Model 1, so any of these models are plausible. 
MRC1349-265: Model 7 was selected by A ln(Evidence) = 
1.9 over Model 1, the highest-evidence model including a 
disk. 

MRC1355-215: Model 7 was selected over the next best 
model, Model 1, by A ln(Evidence) = 3.7. 
MRC1355-236: Model 1 was selected with strong evidence. 
MRC1359-281: There is strong evidence for Model 1. The 
fitted angle is 7° ± }, which does not agree with the angle 
fitted from the first analysis, 4° ± }, for the case in which 
the accretion disk is not fixed at the redshift of the narrow 
lines. 

7.3 Fits to the emission spectra 

Table [T2l gives the results of the model selection for the case 
in which the NLR and the accretion disk are fixed at the 
same redshift. For 14 out of 19 sources, the la error range of 
the best-fit disk angle from the new analysis overlaps with 
the la error range of the best-fit disk angle for the case 
in which the disk is allowed to shift with respect to the 
NLR. The remaining five are MRC0450-221, MRC1114-220, 
MRC1208-277, MRC1222-293 and MRC1359-281. 

The narrow line centres of a number of the models were 
not strongly constrained in the fitting process, and since the 
disk and NLR are fixed at the same redshift in this regime, 
an incorrectly fitted narrow line centre leads to an offset in 
the disk position as well, giving an erroneous fitted angle. 
This potentially affects MRC0327-241 (in which the narrow 
line position differed by 20A between the fit from the analy- 
sis in which the disk shift is allowed and the one in which it 
is not); MRC1114-220, in which there was degeneracy in the 
fitted position of the narrow Ha; MRC0413-296, MRC0430- 
278, MRC1222-293 and MRC1359-281, in which the narrow 



lines were fit with very broad widths (> 30 A), and hence 
their positions may not have been fixed to very good accu- 
racy. This effect might plausibly be the reason that the disk 
angles from the two analyses do not agree for MRC1114- 
220 and MRC1359-281. However, at least MRC0450-221, 
MRC1208-277 and MRC1222-293 have different disk angles 
as a direct result of the fixing of the accretion disk at the 
redshift of the NLR. 

In addition to the inconsistency in disk angles for these 
five sources between the two analyses, the two sources with 
the most extreme changes in fitted disk angle (MRC1208- 
277 and MRC1222-293) have changed from having strong 
Bayesian evidence in favour of a disk, in the regime where 
a disk shift is allowed from the NLR, to having strong 
Bayesian evidence against the presence of disk, when the 
disk is fixed at the redshift of the NLR. The volume of prior 
parameter space was reduced in this second analysis by the 
removal of a variable, which, if this variable was extraneous, 
would have increased the Bayesian evidence for the presence 
of accretion disks. The fact that the evidence remained sim- 
ilar in the majority of cases, but declined sharply in these 
two cases, implies that the variable describing the velocity 
shift between the disk and the NLR is required. 

Figure \W\ shows the fitted disk angles for the best-fit 
models including a disk, in the analysis where the disk is 
fixed at the same redshift as the NLR. Comparing these fit- 
ted angle distributions with those from the case where a disk 
- NLR velocity shift is allowed (shown in Figure , it can 
be seen that these distributions are reasonably similar; only 
MRC1208-277 and MRC1222-293 have fitted angles which 
change dramatically. There is a general tendency for the disk 
angle error bars to be smaller in the new analysis, which is 
due to the reduced number of variables tightening the fit. 
It should be noted that the error bars plotted are simply 
la error bars given the quasar spectrum and the model: no 
error has been added to account for the choice of model. 

7.4 Relationships with radio source size 

Figure [17] shows, for the new analysis, sine disk angle plot- 
ted against projected linear source size, which if the radio 
jets are perpendicular to the disks should be correlated due 
to geometric effects. Fixing the disk at the redshift of the 
NLR has strengthened this from a la correlation, in the case 
that the model allowed a disk velocity shift (see Figure [8]), 
to a 2a correlation (Kendall r coefficient 0.779, probability 
that a correlation exists 97%). This correlation remains at a 
significance level of 2a if the CSS sources are removed from 
the sample (Kendall r coefficient 0.982, probability that a 
correlation exists 96%). 

The source sizes were deprojected using the method in 
Section 16.31 and the resultant sizes are given in Table 1131 
with those from the original analysis for comparison. A his- 
togram of the deprojected sizes is shown in Figure 1181 the 
fitted angles from the new analysis deproject the source sizes 
to a slightly narrower range than the original analysis (see 
Figure [9]) . This does not imply a preference for the mod- 
els with no disk shift, since the range of deprojected sizes 
merely gives an indication as to the typical size to which a 
source expands before it becomes quiescent. 

Figure [19] shows the cumulative distribution of the de- 
projected source size for the analysis of the models with 
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Table 12. Summary of the model fitting results for the analysis in which the accretion disk is fixed at the same redshift as the NLR. 
Columns 1 and 2: Quasar index and MRC name; Column 3: Index of the best-fit model. If the selected model does not include an accretion 
disk, the index of the highest-evidence model which includes an accretion disk is shown in brackets. Column 4 '■ Indication of whether there 
is evidence for a disk according to the odds ratio, with codes as follows: SD = strong evidence for a disk; MD = moderate evidence for a 
disk; WD = weak evidence for a disk; PD = a possible disk, i.e. the best fit is a non-disk model, but there is only inconclusive, weak or 
moderate evidence for this; ND = strong evidence against the presence of a disk; Columns 5 and 6: Natural logarithmic evidence difference 
between the best-fit model, in cases where this includes an accretion disk, and the second best model (Column 5, Aln£ 2 d ), and in cases 
where the second best fit also includes a disk, between the best- fit model and the best- fit model without a disk (Column 6, Alni?^); 
Columns 7 and 8: Natural logarithmic evidence difference between the best- fit model, in cases where this does not include an accretion 
disk, and the second best model (Column 7, Aln Ef n ), and in cases where the second best fit also does not include a disk, between the 
best-fit model and the best-fit model with a disk (Column 8, Alni?^); Column 9: Best-fit angle of the disk and error in degrees. For 
analysis in which the disk is allowed to shift with respect to the NLR — Column 10: Index of the best-fit model; Column 11: 
Indication of whether there is evidence for a disk according to the odds ratio, codes as before; Column 12: Best-fit angle of the disk and 
error in degrees. 



the disk fixed at the redshift of the NLR. This distribution 
is marginally more consistent with the comparison models 
of uniformly expanding sources than the distribution found 
in the original analysis (see Figure llO]) . Both distributions 
are consistent, within the limits of small number statistics, 
with the picture of sources expanding uniformly up to some 
cut-off size, though there are indications in both cases of an 
excess of small sources. 

7.5 Relationships with radio luminosity 

The sine best-fit disk angles for the analysis with the ac- 
cretion disks fixed at the NLR redshift are plotted against 
the 178 MHz radio luminosity in Figure 1201 These angles 
are consistent with the envelope of angles predicted by the 
receding torus model. However, the correlation seen be- 
tween the sine of the disk angle and the low-frequency ra- 
dio luminosity in the original analysis, at a significance of 
> 1.5cr according to Kendall's r test, has now fallen to 
< 0.5a (Kendall's r coefficient 0.105, probability of corre- 
lation 25%). This illustrates the strong dependence of the 
analysis results on the model assumptions, such that a large 
change in fitted angle for only two quasars is enough to dis- 
rupt this correlation. 

In order the test the fragility of the luminosity - disk 



angle correlation, the statistical analysis was repeated 24 
times, with the disk angle modified in each case by a Gaus- 
sian scatter with standard deviation of 2.16, the standard 
deviation of the distribution of differences between the disk 
angles with the disk fixed at the redshift of the NLR, and 
the disk angles fitted when the accretion disk was allowed a 
velocity shift. Only 8 out of 24 test correlations have higher 
significance than the original correlation between luminosity 
and disk angle of a disk not fixed with respect to the NLR, 
which has a significance of 1.8a; however, 13 out of 24 test 
correlations have a significance of > 1.5a. A larger sample 
of quasars might go some way towards confirming whether 
this is a real correlation, but it is clearly necessary to fully 
test a wider range of emission models. 



7.6 Distribution of angles 

The cumulative probability distribution for the disk angles 
and the theoretical distribution of jet angle for different val- 
ues of the Lorentz factor is shown in Figure 1211 With the 
assumption that the axis of the accretion disk is coincident 
with the radio jets, then the distribution of disk angles is 
consistent with a value of the jet Lorentz factor around 
7 ~ 20, as for the previous analysis (see Section fo.5p . 
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Figure 16. Disk angles for model fits with an accretion disk, for 
the analysis where the disk emission is fixed with respect to the 
narrow lines. The black points indicate the best-fit disk model for 
each quasar. Circular symbols indicate evidence for a disk, square 
symbols indicate a possible disk (inconclusive evidence, or weak 
to moderate evidence against a disk), while triangular symbols 
indicate strong evidence against the presence of a disk. White 
points mark fits which fall within one Jeffreys' criterion of the 
best fit, and these can be seen in all cases to be consistent with 
the best-fit disk angle. All angles derived from non-best-fit disk 
cases are plotted in pale grey. 



7.7 Discussion 

The fitted disk angles have proved to be only partially robust 
to a change in the model in which the accretion disk is fixed 
at the same redshift as the NLR, although the majority were 
consistent. The small number of objects in this sample make 
it difficult to judge whether the disappearance of the disk 
angle - low-frequency radio luminosity correlation in this 
new analysis occurs because the initial results are a fluke, or 
is a genuine effect from the alteration of the models. Since it 
is not clear that the narrow lines and accretion disk should 
be fixed at the same redshift, this new analysis should serve 
as a measure of the fragility of the results with respect to a 
change in the model. 



8 CONCLUSIONS 

All bar one of the 19 quasars are best fitted with an emis- 
sion model which includes more than one component of 
broad-line emission, in both analyses; there is overwhelm- 
ing evidence that in this sub-sample of quasars, a simple 
one-component broad-line emission model is insufficient to 
describe the physical processes at work. There is strong 
Bayesian evidence that ten out of nineteen quasars from the 
analysis in which a shift is permitted between the BLR and 
NLR possess accretion disks, and all but one are consistent 
with the emission model of a circular Keplerian accretion 
disk in addition to a single-peaked, symmetric emission line 
arising from a separate BLR. The one quasar for which this 
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Figure 17. Sine best-fit disk angles versus the projected lin- 
ear size of the quasar in kpc for the analysis in which the disk is 
fixed at the redshift of the NLR. The disk angle error bars include 
both the errors returned from the individual fits and a systematic 
error calculated from the variation in results from using differ- 
ent random seeds to initialise the Monte Carlo engines, added in 
quadrature, but do not take into account any error due to model 
selection. The non-CSS FR1I sources are plotted with circles, and 
CSS sources are plotted with squares. The two sources plotted 
with crosses are core-dominated sources without measured pro- 
jected sizes. The black points have weak to strong evidence for a 
disk, grey points have possible disks, and the white points have 
Bayesian evidence against the presence of a disk. Index numbers 
are matched to source names in Tabled The black line shows the 
best linear fit to all sources excluding the core- dominated sources. 
There is a 2a correlation by Kendall's r test between the sine of 
the disk angle and the projected size both for all sources exclud- 
ing the core-dominated sources, and for the 11 non-CSS FRII 
sources. The dotted, dot-dot-dot-dashed, dot-dashed, dashed and 
long-dashed lines show loci of constant deprojected size, for 10 
kpc, 25 kpc, 100 kpc, 1 Mpc and 5 Mpc respectively. 



does not apply appears to have very complex broad emis- 
sion, consisting of at least three components. If the model 
does not allow for a shift between the BLR and NLR, strong 
evidence for disks is seen in eight out of the nineteen sources, 
and all but three are consistent with the presence of a cir- 

cula r accretion disk. 

lEracleo us & Halpcrn dl994h . lEracleous fc Halpernl 
[j2003r i and Istrateva et all l|2003n found that only 3 - 20% 
of AGN have obvious double-peaked lines (see Section E^}, 
in contrast to the ~ 84 - 95% of this sample of quasars 
which were found to be consistent with the presence of a 
thin, optically-emitting accretion disk, despite only a few of 
the line profiles appearing obviously double-peaked to the 
eye. It may be that accretion disk emission has been seen 
in so few AGN because single-peaked broad emission lines, 
arising from fast-moving clouds located outside the thin 
disk or from a wind originating from the disk itself, have 
a tendency to swamp the disk emission. In these cases, a 
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log( Source size deprojected with disk angle / kpc ) 

Figure 18. The solid black line shows source sizes deprojected 
with the fitted disk angles from the models where the disk emis- 
sion is fixed at the same redshift as the NLR (see Table I13jl . 
The source sizes from deprojecting with the fitted angles from 
the analysis in which the disk emission is allowed to shift with 
respect to the NLR are shaded in grey. The projected sizes are 
plotted with a dotted line. Note that the two core- dominated 
sources have no measured size and are therefore excluded. 
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Figure 19. Cumulative distribution of source sizes deprojected 
using the best- fit disk angles from the models where the disk emis- 
sion is fixed at the same redshift as the NLR (solid line). The 
projected source sizes (dotted line) and the source sizes depro- 
jected with the angles from the analysis where the disk emission 
is allowed to shift with respect to the NLR (grey line) are shown 
for comparison. The curves show the predicted distribution if the 
heads of the sources are expanding at a constant rate; the dot- 
dashed line is normalised at the deprojected size of the largest 
source in the sample, the dashed line is normalised at 1 Mpc. 
Note that the two core-dominated sources have no measured size, 
and are boosted into this sample by strong core emission, and are 
therefore excluded from this plot. 
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Table 13. Summary of projected and deprojected source sizes. 
Columns 1 and 2: Quasar index and MRC name; Column 3: 
Source type — CD indicates sources with core-to-lobe radio flux 
ratio at 10 GHz of greater than 1, CSS indicates sources with 
Qopt > 0.5 and projected source sizes < 25 kpc, FRII indicates 
non-CSS FRII, which encompasses all other sources; Column 4- 
Projected source sizes in kpc (sources of these values are given in 
Table |SJ; Column 5: Deprojected source sizes in kpc, calculated 
using the fitted disk angles from the models in which the accre- 
tion disk is fixed at the same redshift as the NLR; Column 6: 
Deprojected source sizes in kpc, calculated using the fitted disk 
angles from the models in which the accretion disk is allowed to 
have a different redshift from the NLR. 



complex model and high signal-to-noise ratio spectra are 
required in order to retrieve the accretion disk contribution 
to the emission. It is therefore essential to extend this 
analysis to a larger sample of low-radio-frequency-selected 
quasars; optically-selected quasars are predicted to cover 
a narrower range in disk angles, as even lightly reddened 
quasars (at ~ 45°) are much more likely to fall through the 
magnitude limit of a survey, and the optical emission may 
be angle dependent. 

There is a possible correlation of the fitted disk angles 
with the projected source size, in agreement with the ex- 
pected projection effects if the accretion disk is perpendicu- 
lar to the jet. There are two CSS sources which do not appear 
to lie on the same relation as the other sources: MRC0222- 
224 seems to be an intrinsically smaller and more reddened 
source than the others, while MRC 11 14-220 may either be 
an outlier, or have a precessing or misaligned disk. When 
deprojected using the best-fit disk angles, the source sizes 
are consistent with a model in which the accretion disk is 
perpendicular to the radio jets, and the heads of the radio 
jets expand at approximately constant speed up to a size 
~ 1 Mpc; the paucity of sources larger than this limit can 
be explained by the average duty cycle of these AGN. 

The distribution of fitted disk angles is consistent with 
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• / ■ : Non— C55 / CSS sources with evidence for disks 
~ O / □ : Non-CSS / CSS sources with possible disks 

O / □ : Non-CSS / CSS sources without evidence for disks 
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Figure 20. Sine best-fit disk angles versus logarithmic 178 MHz 
luminosity. The disk angle error bars include both the errors re- 
turned from the individual fits and a systematic error calculated 
from the variation in results from using different random seeds 
to initialise the Monte Carlo engines, added in quadrature, but 
do not take into account any error due to model selection. The 
non-CSS FRII sources and core-dominated sources are plotted as 
circles, and CSS sources are plotted as squares. The black points 
have weak to strong evidence for a disk, the grey points have 
possible disks (there is inconclusive evidence, or weak to mod- 
erate evidence against a disk) and the white points have strong 
evidence against a disk. Index numbers are matched to source 
names in Table [T] The best-fit line (minimising \ 2 ) is shown in 
black. The dashed line sho ws the r a dio luminosity-dependent crit- 
ical angle ( following [Willott et all l|2000h l with fiducial angle 45° 
at logLi7gMHz= 27, and minimum quasar fraction 10%. 



the calculation of expected jet angles for Doppler-boosted 
jets with Lorentz factors of ~ 20, if the accretion disk rota- 
tion axis and radio jets are assumed to be coincident. The 
calculated jet angle distribution may match the disk angle 
distribution better if a more physically reasonable model was 
used for the jet, with a core of higher Lorentz factor than 
the surrounding jet sheath, as opposed to this simple single 
Lorentz factor model. 

There is a weak correlation (> 1.5a significance by 
Kendall's r test) between the low-frequency radio luminos- 
ity and the sine of the angle between the accretion disk axis 
and the l ine of sigh t . This is predicted by the receding torus 
model of lLawren cc (1991) if the accretion disks are perpen- 
dicular to the radio jets, such that the disk angles match 
the jet angles. In this model, the dust is sublimated by the 
AGN nucleus, leading to larger opening angles of the obscur- 
ing dusty torus for more radio-luminous sources, and hence 
quasars are seen up to larger jet angles. The largest accre- 
tion disk angle measured is 48° , which is consistent with the 
opening angle predicted for powerful FRII radio sources. A 
larger sample size is needed to confirm the radio luminosity 
- disk angle correlation, but this is vital, as it would yield 
the first direct test of the receding torus model. It should be 
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Figure 21. The cumulative distribution of the best-fit disk angles 
from the analysis in which the disk redshift is fixed at the redshift 
of the NLR (solid black line). The grey line shows the cumulative 
distribution for the best-fit disk angles for the analysis in which 
the accretion disk is allowed to shift with respect to the NLR. The 
curves s how the theoreti cal distributions of jet angle, modelled 
from the I Jackson fc Walll l|l999l l relation and normalised up to an 
angle of 45°, for different values of the Lorentz factor: the dashed 
line shows the expected cumulative distribution for 7 = 3, the 
dot-dashed line for 7 = 8.5, and the dotted line for 7 = 20. The 
vertical dot-dot-dot-dashed line indicates a fitted disk angle of 
45°. 



noted that when the analysis was performed again with the 
models in which the accretion disk was fixed at the redshift 
of the NLR, this result disappeared, and so the radio lumi- 
nosity - disk angle correlation is strongly dependent on the 
model used. 

A study of the velocity shifts between the line compo- 
nents revealed the single-peaked broad lines to be formed 
at similar redshifts to the narrow lines. The disk emission 
tended to be redshifted with respect to these lines for high- z 
sources and blueshifted for lower- 2 sources. Since the accre- 
tion disk is expected to be at the systemic velocity, this 
would imply that the single-peaked broad lines and the nar- 
row lines are formed in outflows for more luminous, high- 
z sources, and in infalling clouds for closer, less luminous 
sources. Such a scenario could arise if powerful winds were 
causing outflows of gas from the high-z, luminous sources, 
while these winds have stalled in the less-luminous quasars 
seen at z < 1, resulting in infall of the gas back towards the 
galaxy nucleus. 

The results of this paper are all dependent on the va- 
lidity of the disk model used in the Baye sian fitting. This 
model, taken from [Ch en fc Halpernl (|l989h . describes a cir- 
cular, Keplerian disk. The illuminating flux from the central 
black hole was modelled as falling off with a radial exponent 
of 3. In order to investigate the effects of range of permuta- 
tions to this basic model, such as a variety of disk emissivities 
and warped or elliptical disks, greater spectral resolution, 
a larger quasar sample size, and in particular, multi-epoch 
data to study the variability of the emission line profiles are 
required. 
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In addition, similar tests were made in which the 
Bayesian fitting parameters were unchanged, but the spectra 
were altered. The first of these tests added Gaussian random 
noise to the spectrum, with FWHM of 5% of the data values, 
in order to see qualitatively what effect greater noise has on 
the fitting process. For MRC0450-221, the posterior prob- 
ability for the local velocity dispersion occupied a similar 
region in parameter space to the total space occupied by all 
fits in the test with different Bayesian random seeds; it ap- 
pears that the velocity dispersion parameter is constrained 
to a fairly large region of parameter space, and within this, 
for some tests, local minima are found. For MRC0437-241, 
the parameters constrain in similar regions independent of 
the noise added. 

The second of these tests was the fitting of spectra 
smoothed with boxcar widths of 3, 5, 7 and 9 pixels, where 
each pixel has a scale of approximately 2 A. For MRC0437- 
244, the smoothed spectrum allowed the probability to con- 
strain marginally more tightly, but in the same regions. This 
was also true of MRC0450-221 for the boxcar smoothing 
widths of 3, 5 and 7 pixels. At a boxcar width of 9 pixels, 
the posterior probability distributions of the fit alter signif- 
icantly, as real structure has been lost from the spectrum. 



APPENDIX A: TESTS FOR STABILITY OF 
THE FITTING 

Detailed tests were run on two of the quasars selected at ran- 
dom (MRC0437-244 and MRC0450-221), in order to check 
the Bayesian routines for stability. In each set of tests, one 
of the BayeSysS parameters was changed: the Bayesian ran- 
dom seed; the rate at which the Monte Carlo exploration 
focuses in on the region of maximum likelihood; and the en- 
semble of parallel explorations of the parameter space during 
the Monte Carlo optimisation. Plots of the posterior proba- 
bilities of the disk parameters were created, each test being 
overlaid in a different shade, to identify fits in which the 
posterior probability converged in a different region to the 
majority of the fits. 

An example of one of the stability tests is shown in Fig- 
ure [4] this shows the convergence of some of the disk param- 
eters for MRC0450-221, as well as the correlations between 
these parameters, when the Bayesian fitting is initialised 
with five different random seeds. The posterior probability 
distribution for each seed is plotted in a different shade of 
grey. The local velocity dispersion, Av, of the disk mate- 
rial converges in two slightly different regions of parameter 
space, depending on the Bayesian seed, though these areas 
are within a reasonable error bound of each other. The local 
minima do not have vastly different evidence values in this 
case (A In (Evidence) ~ 2.5 between the highest and lowest 
evidence cases); this is due to the fact that the profile of 
the disk e mission varies only we akly with the local velocity 
dispersion l)Strateva et al.l [20031 ). The variation in velocity 
dispersion was similar through all of the parameters tested: 
all other parameters were shown to be stable in each of the 
tests for MRC0450-221. All parameters were stable in the 
tests for MRC0437-244. 
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